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I. INTRODUCTION 



In classical mechanics, determining the motion of three bodies which attract each other 
according to Newton's universal law of gravitation is the most celebrated of all dynamical 
problems (see, for example, jH, 0). Over the period since 1750, it has attracted the attention 
of some of the greatest mathematicians and, in particular, Euler, Lagrange, Jacobi, Poincare 
and Levi-Civita. It cannot be solved exactly. However, recently the advent of high-speed 
computers has opened up a whole new approach to the problem by making possible step- 
by-step numerical integration of the differential equations of motion from the initial time 
to any desired later time. The quantum three-body problem also has a rather well-known 
history, in particular for systems governed by Coulomb forces. Some contributions will be 
reviewed in this article. 

Suppose we have a quantum system made up of three particles, two having the same 
charge and the third having a charge of the same magnitude but of opposite sign, interacting 
only through Coulombic forces. According to non-relativistic quantum mechanics, such a 
system may, or may not, form bound states, depending on the relative masses of the particles. 
For instance H^(p, p,e~), Ps~(e"*", e~, e~) or H~(p,e~,e~) are stable, while (p, e"'",e~) or 
(p, p, e~) are not. The critical values of the relative masses at which bound states appear 
have not been established. However, we shall describe in this article some of the interesting 
and often ingenious methods which have been devised to obtain bounds on the number of 
bound states of a given system and also on the critical mass ratios at which such states 
appear. Our aim is also to describe the properties of the stability domain as a function of 
the masses or inverse masses of the particles involved. 

We then extend our study in several ways. First we consider three-body systems with 
arbitrary charges. There are sets of charges for which stability is obvious, i.e., is achieved 
whatever values are assumed for the masses. On the other hand, there are more interesting 
situations, where stability requires very large values of certain mass ratios. 

Another field deals with systems containing more than three unit charges. In the case 
of four charges, whose best known prototype is the hydrogen molecule (p,p,e~,e~), we 
discuss for which masses there exists at least one bound state stable against spontaneous 
dissociation. 

The sections which follow will be devoted to systems with five, six or more unit charges. 
Here, the boson or fermion character has to be specified when dealing with more than two 
identical particles. 

We will also repeat some of the studies with only d = 2 space dimensions. In the one- 
body case, (or for two self- interacting particles), the situation is dramatically different for 
d < 2 and d > 2 since in the former case, any attractive potential, however weak, will 
produce a bound state. One expects also qualitative changes for more complex systems. 
Some conclusions, and a list of unsolved problems will be included in the last section. 

Some useful results on Schrodinger operators, an introduction to methods for solving the 
few-body problems, in particular to the Born-Oppenheimer approximation and variational 
calculations, and a recollection of elementary calculations used historically to establish the 
stability of 3- and 4-body systems, are given in Appendixes. 

This review contains both a summary of rigorous results and a survey of accurate numer- 
ical methods. The knowledge of precise binding energies of few-charge systems is extremely 
useful in comparing the efficiency and accuracy of various quantum mechanical methods. 
It also helps development of new methods and the improvement of those already in exis- 



tence, as well as the understanding of electron-electron, electron-positron or electron-lepton 
correlations and dynamics. 

The investigation of the stability of few-charge systems is not only a very good test- 
ing ground for different models and methods but it also leads to important applications. 
The study of stability provides very useful information for various physical systems, e.g., 
positronic atoms 0,0iSB0], antimatter compounds systems of atoms and antiatoms 
0, E^, or charged excitons (system of electrons and holes) in semiconductors [lH. New 
experimental investigations of the fine structure of charged excitons in semiconductor 
quantum dots have shown the existence of bound states of three electrons and a positive 
charge or four electrons plus two-positive charges. The stability studies presented in this pa- 
per clearly show that such systems are not stable if only the pairwise Coulomb (or Coulomb 
like) interaction is included and their existence is therefore due to some other binding mech- 
anism, e.g., confinement. There is also an interesting analogy between multiquark systems 
in hadron physics, and small exotic molecules: the property of flavour independence states 
that the same potential applies to quarks of different masses, in the same way that the same 
Coulomb interaction is felt by electrons, muons or antiprotons. For such systems, similar 
patterns are observed for configurations leading to maximal stability against spontaneous 
dissociation 

In the preparation of this review article, it has been necessary to read hundreds of papers. 
Unfortunately, it is impossible to discuss or even list all the contributions that have been 
made to the subject of few-charge systems in quantum mechanics. For example, there 
are more than twenty papers dealing with the binding energy and static properties of the 
positronium ion, with considerable overlap between them. In each section, the key pioneering 
papers are mentioned as well as a wide range of more recent contributions which should make 
it possible for the reader to explore the parts of the literature that are not referred to. We 
apologise in advance to any author of a significant contribution that has inadvertently been 
omitted. 



II. STABILITY CRITERIA 



As we are considering self-interacting, translation-invariant Hamiltonians in this review, 
we shall assume that the centre-of-mass motion is removed (see Appendix El, leading to the 
"interaction Hamiltonian" -f^int- In the case of three unit charges [mf , m2 , m^) , one can 
assume m2 < m^. Then, in simplified units, 

= -i,,, = -i^^ . (2^1) 

nil + fT^3 

the ground-state energy of the hydrogen atom-like system made up of particles 1 and 3, can 
be shown to be the lowest continuum threshold of ifint- For systems with more than three 
charges, there is always such a threshold at which the system separates into two or more 
smaller neutral or charged systems made up of one or more particles. 

Any eigenfunction with eriergy below this threshold will be square-integrable and thus 
correspond to a bound state |l3| . 

This behaviour will be called genuine stability or simply, stability. Of course, it is sufficient 
to establish that some upper bound, as given for instance by the Rayleigh-Ritz variational 
method (see Appendix^, lies below i?th- However, if it cannot be established using the 
Rayleigh-Ritz method that an upper bound exists that is below i^th, this does not necessarily 
mean that the state under consideration is not bound. It may be that the nature of the 
variational trial function makes it impossible to obtain a sufficiently good upper bound. 
One possible way of showing that a state is not bound, when the Rayleigh-Ritz method is 
inconclusive, is to obtain a lower bound to the energy of the state and show that this is 
above Eth (see Appendix \K\ Section 1X3]) . 

Except in a few cases, we shall not consider the possibility of bound states in the contin- 
uum. The existence of such states is conceivable but they are rather pathological and are 
expected to be unstable under perturbations fl^. 

Forms of stability weaker than genuine stability can also occur. 

1. Imagine two heavy nuclei or particles surrounded by electrons. The particles behave 
almost classically in the Born-Oppenheimer potential V{r). If V{r) has a local minimum at 
internuclear distance r = tq, classical particles with sufficiently small kinetic energy would 
remain in the vicinity of this minimum, even if its value V{rQ) lies above the dissociation 
energy. We shall call this behaviour classical stability. 

2. Imagine a 3-body state (a^,6~,c~) with angular momentum and parity = 
this not including intrinsic spins and parities. Then the decay 

(a+,r,c-) ^ (a+,r)i, + c- , (2.2) 

cannot take place if the atom in the final state is in its Is ground-state with angular mo- 
mentum and parity = O"*". The decay requires higher states, for instance 2p. This means 
that the threshold is higher, as long as radiative corrections are neglected. Otherwise, the 
reaction 

(a+,r,c-) ^ (a+,r)is + c- + 7 , (2.3) 
would have to be considered. We shall call this behaviour metastability. 



III. PARTICULAR THREE-UNIT-CHARGE SYSTEMS 



In this section, we present the methods apphed to examine the stabihty of three unit- 
charge systems, and discuss in some details configurations which have been extensively 
studied. The next sections will be devoted to global studies of the stability region in the 
parameter space. 



A. Notation 

We consider systems with one charge, qi, of mass mi, and two charges opposite to qi, 
say q2 = qs = —qi, with masses 1712 and ms, respectively. By charge conjugation invariance, 
it does not matter whether gi > and q2 = q-s < or vice-versa. By scaling, one can take 

qi = -q2 = -qs = i- 

The Hamiltonian thus reads 

^=^ + ^ + ^ + ^-^-^, (3.1) 
2mi 2m2 2m3 ri2 

where r^j = \rj — ri] is the distance between particles i and j. 

If we separate off the centre-of-mass motion, we obtain the interaction Hamiltonian 

i^in. = H- ^P^^P^ + P^^\ . (3.2) 
2(mi + m2 + ms) 

Then if we take the origin of -f^int to be at particle 3, 

^int = -^V^^ - -^V^, - —v., . v., - 1 + 1 - — , (3.3) 
2/ii 2fi2 ms n r2 ru 

where ri and r2 are the position vectors of particles 1 and 2, respectively, with respect to 
the third particle as origin, m^ is the mass of the ith particle, and 

mima m2m3 

/^i = ] , = ■ , 3.4 

mi + ms m2 + ms 

where fii, for example, is the reduced mass of particles 1 and 3. 

Hint represents the motion of two particles (more properly, pseudo-particles) of mass /ii 
and /i2, respectively. It differs from the usual form for the Hamiltonian for two particles 
as it contains a coupling term between the momenta of the particles. This is known as the 
Hughes-Eckart or mass-polarisation term jisl . iH 17 . 



If ms is very much larger than mi and m2, the effects introduced by locating the origin 
at particle 3 are very small. For this reason, they are normally neglected in the analogous 
many-body problem of an n-electron atom. 

The mass-polarisation term can be eliminated by using Jacobi coordinates, first intro- 
duced by Jacobi and Lagrange, in the context of the three-body problem in Newtonian 
mechanics 0. If we transform to new internal coordinates, (see Appendix IX]). 



which corresponds to taking the origin of the coordinates of particle 2 to be at the centre of 
mass of particles 1 and 3, H^^t is of the form 3, 1^ • 



gint = ^V^--^V^+ \ --- ^ (3.6) 

where 

m2 (mi +7723) ms , mi 

/^x = , ^iy = 7 . ■. r , V = ■ , 7] = ■ ■ 3.7 

(mi + m2 + m^) mi + m^ mi + m^ 

Note that though x = ri, differs from as y, and not r2, is held constant during the 
partial differentiation. As seen in Appendix^ using Jacobi coordinates simplifies the kinetic 
energy part of ifint at the expense of making the potential energy part more complicated. 



B. H+ 

The hydrogen molecular ion, H^, belongs to a special category of three-body systems that 
contain two identical particles and a third particle with charge of equal but opposite sign. 
Using the comparison theorem, see Appendix Hill jl^] was able to prove that systems 
of this type always have at least one bound state. He was able to show further that these 
systems have only one bound state if mi/ ma < 0.210101636, where m^ is the mass of the 
identical particles and is the mass of the remaining particle. 

In the case of H^, where the two identical particles are protons and the remaining particle 
an electron, mj/m^ = 1836.2. Thus this system is far outside the region in which it can be 
shown that only one bound state exists. As might be expected, it can be shown that this 
system has many bound states. 

The large difference in mass between the protons and the electron means that the Born- 
Oppenheimer approximation (Appendix [H]) gives accurate results when applied to this sys- 
tem. In this approximation, the nuclei are fixed and the energy of the electron is calculated 
as a function of the internuclear distance, R. The energy of the system is then obtained by 
calculating the energy of the nuclei moving in a potential made up of the electronic energy 
and the internuclear potential. 

The Hamiltonian for the electron moving in the field of the fixed nuclei is 

^ = --^V^----, (3.8) 

2me ta tb 

where m^. is the mass of the electron and ta and tb are the distance of the electron from 
nucleus A and 5, respectively. The units are such that h = 1 and e = 1, where e is the 
charge on the proton. In terms of prolate spheroidal (confocal elliptical) coordinates A, fi, 
(p (see, for example, jl^]), the potential, V, in this Hamiltonian is of the form 

where 

rA + TB ta-tb 
— R — ' — R — ' 



and is the azimuthal angle with respect to Cartesian coordinates with origin at the mid- 
point of the internuclear axis and 2;-axis directed along this axis. 

As is well known (2^ |23| , the Schrodinger equation with this Hamiltonian is separable 
in terms of prolate spheroidal coordinates. Details of papers describing how the resulting 
ordinary differential equations can be solved to give physically acceptable solutions are listed 
on page 2 of Ref . [2l| . Details of square-integrable wave functions and energies are given by 
Bates et al. jl^. More accurate results for the ground-state energy have been obtained by 
Wind I25I] and Peek 0|. For a review of calculations of this type, see Bates and Reid [27^ . 

Calculations of bound-state energies for Hg" and the very closely related system HD+, 
that include corrections to the Born-Oppenheimer approximation, have been carried out 
by Beckel et al. [HI, Wolniewicz and Poll liol. Wolniewicz and Orlikowski js^]. Moss jsH, 



Taylor et al. Korobov [33], Hilico et al. [34] and Frolov J35]. The results of Wolniewicz 
and Poll's calculations for HD+ are compared with experiment by Carrington and Kennedy 
[2^, who also include a detailed discussion of the underlying theory and an extensive list 
of references to earlier work on HD+. Theory and experiment are in agreement that there 
are a large number of bound states for these systems. The bound states are classified using 
a vibrational quantum number, f , and a rotational quantum number, J. See, for example, 
Bransden and Joachain v may have a value as high as f = 19 for and v = 21 for 
HD"*". In general, each value of v is associated with a number of values of J, starting with 
J = 0. 

Recently, Carbonell et al. found a new state of H^, with J = 0, but with an antisym- 
metric coordinate-space wave-function for the two protons. The binding energy is extremely 
small, of the order of 10~^a.u. This state has been discovered by studying how the a"'"H 
scattering length evolves as a function of the mass of the projectile a"*". The case of 
corresponds to a being a proton. The scattering equations at zero energy are solved in the 
configuration-space Faddeev formalism. The energy of the new state of was confirmed 
by a very accurate variational calculation. 

(p,e~,/i"'") and (/i"*", e~, /i"*") are systems similar to but with one or both protons re- 
placed by a positively charged muon. The mass of the muon is 206.8 me. This is much less 
than the mass of the proton, but sufficiently greater than mg that the Born-Oppenheimer 
approximation gives accurate energy values when applied to these systems. In this approx- 
imation, the internal energy of (p,e~,/x"'") or e~, /i"*") is obtained by calculating the 
energy levels of a pseudo-particle with mass equal to the reduced mass of (p, or (/i"*", /^+) 
moving in the same electronic and Coulombic internuclear potential as in the case of 
and HD"*" (see Appendix IH]) . The reduced masses of {p,fi'^) and (/i"*", /x"*") are 185.9 mg and 
103.4 me, respectively. These values are much less than 918.1 mg, the reduced mass of (p,p). 
It follows from this that these systems do have bound states but substantially fewer than 
in the case of Hj. See, for example, Bransden and Joachain 17]. For a non-adiabatic cal- 
culation of f/x"*", e~, M"*"), see, e.g., j^^, and for the asymmetric (p, e~, /x"*"), (d, e",/!"*") and 
(t,e~,/x"'") systems, 



C. H- 

H" contains two electrons and a proton and is thus another example of a three-body 
system that contains two identical particles of unit charge and a third particle with equal 
but opposite charge. In this case, mj/m^ = 0.00054. It follows from the result of Hill \2(\ 
referred to earlier that only one bound state of this system exists. This is the singlet ground 



state. 

It has been known from early on in the development of quantum mechanics that the singlet 
ground state is bound. See Appendix |D| Variational upper bounds of sufficient accuracy 
to demonstrate this were obtained by Bethe Hylleraas 0| and Chandrasekhar j4S 
Very accurate upper bounds for the energy of this state have been obtained by Pekeris 
Frolov and Yeremin [i^l , Cox et al. ji^ and Drake et al. ji^ . 

There exist a number of autodetaching or autoionising states of H~ below the continuum 
thresholds corresponding to excited states of the H atom. See Massey [47], Bhatia and 
Temkin ^] and Bhatia |4^. In particular, there is a triplet P° state below the n = 2 
continuum threshold. This is not a bound state as it is embedded in the n = 1 continuum; it 
has a small but non-zero width indicating that autoionisation can occur through interaction 
with states in this continuum. 

__There is also a triplet F** state with energy 0.0095 eV below this threshold. See Bhatia 
and Drake jHH . In this case autoionisation is forbidden if spin-orbit coupling is neglected 
as in this review article. This is because the n = 1 continuum does not contain any states of 
the required angular momentum and parity for this process to occur. This triplet F'^ state 
is often described as a bound state. See, for example, Drake [5l[ and Ho 152:]. However, 
it should more properly be called a metastable state as it does not satisfy the criterion for 
stability specified in Sec. |n]of this review, namely that its energy must be below the energy 
of the lowest continuum threshold of the system. This is also the criterion used by Hill j20| . 
There is some confusion over this point in the review article on the stability of three-body 
atomic and molecular ions by Armour and Byers Brown ^]. 

Mu~ is a similar system containing a positively charged muon and two electrons. The 
mass of the muon is 206.8 mg. Thus in this case rrii/md = 0.0048 and the system, like H~, 
can only have one bound state. Mu~ has been observed by, for example, Kuang et al. fs^. 
Its ground-state energy has been calculated by, for example, Bhatia and Drachman [551] , 
Fetelenz and Smith jsy] and Frolov and Yeremin [ST^ . 

For a review on two-electron atoms and ions, including H~ and its role in astrophysics, 
see Hi, m, 11" 



D. Ps" 

Fs~, i.e., (e~, e"*", e~), is a system that contains two identical particles, in this case elec- 
trons, and a third particle, a positron of the same mass and equal but opposite charge. It is 
thus the special case of the systems considered earlier for which rrii/md = 1. It follows from 
the result of Hill that this system must have at least one bound state. The existence 
of a bound state of this system was predicted by Wheeler and shown to be the case by 
Hylleraas It was first observed experimentally by Mills [62]. 

The positron is the antiparticle corresponding to the electron. As first shown by Dirac 
jHsj ). it can annihilate with an electron to form gamma rays. The rate of this annihilation 
process for Fs~ was first measured by Mills (H^ ]. 

Details of calculations carried out on Fs~ are given by Ho [s^l and Cox et al. jH^]. See, 
also, Fatil [H^. Very accurate calculations have been carried out for the ground state of 
Fs- by Bhatia and Drachman 0, Ho 0,0, Frolov El Izi| , Fetelenz and Smith 0, 
Frolov and Yeremin 0], Cox et al. jH^, Krivec et al. [T^ and Drake et al. [46*]. 

All the evidence indicates that Fs~ has only one bound state. It is of interest to note that 
all attempts to obtain an energy for a triplet F'' metastable state below the n = 2 threshold 



for Ps^ have been unsuccessful. See Mills j73[ and Ho |52| |. 

Bhatia and Drachman jo^l have examined the behaviour of the energy of the triplet 
state as m^/mj is varied. They show that this energy is below the energy of the n = 2 
continuum if md/rrii > 16.1 or md/rrii < 0.4047. This agrees with the findings described 
earlier for H~ and shows that the energy of this state of Mu~ must be below the energy of 
the n = 2 continuum. Note that rrid/mi = 1 for Ps~. 



E. (p,e-,e+) 

(p,e^,e^), an H atom and a positron, is an interesting case for which the variational 
method fails. Inokuti et al. |7^] were able to show by this method that if the mass of the 
proton is taken to be infinite, a "positron" with mass > 7.8 me would form a bound state 
with an H atom. Frost et al. [t^ reduced this upper bound to 2.625 me by using a more 
flexible trial function. It was reduced still further to 2.20 mg by Rotenberg and Stein f76|. 
They made their trial function even more flexible by including basis functions suitable for 
representing a weakly bound positron in a potential, V{r), with the appropriate asymptotic 
form for an H atom and a positron [itI . IttI I , 

^('■) -|i • (3,11) 

where a is the dipole polarisability of the H atom and r is the distance of the positron from 
the proton. 

The method described at the end of the section in Appendix^on the comparison theorem 
was used by Spruch [t^ to devise a method for showing that no bound state of the e"'"H 
system exists if the mass of the proton is taken to be infinite. 

If we take particle 1 to be the positron and particle 2 to be the electron, Hi^t can be 
expressed in the form 

^int = -TT^V^^ + - + He, (3.12) 

2mi ^ ri 

where 

= -T^'^l, - - - — , (3.13) 
2m2 r2 ri2 

is the Hamiltonian for an electron in the field of two equal, fixed positive charges, one at 
the origin and the other at Vi. As discussed in the Section above on H^, is just the 
Hamiltonian which determines the electronic potential energy of the H2' ion in the Born- 
Oppenheimer approximation. The associated Schrodinger equation is separable in prolate 
spheroidal coordinates. The eigenvalues of He are functions only of ri and the separability 



of the Schrodinger equation makes it possible to calculate them to high accuracy [2J, |2 
Let us consider the adiabatic Hamiltonian 

^ad = i{r2)Hp , (3.14) 

where 

H^ = -^Vl^ + - + Eo{n), (3.15) 
2mi ri 

/(r2) is the unit operator for allowed square-integrable functions of r2 and -Eo(^i) is the 
ground-state eigenvalue of Hg- It is referred to as adiabatic as the potential term, £'o(ri). 



is calculated by fixing ri, then calculating the ground-state energy of the electron in the 
resulting two-centre attractive Coulombic potential. 

It is easy to show that Hi = and H2 = -ffint, where ^^ad and -f^int are as given above, 
satisfy the required conditions for the comparison theorem and have the same continuum 
threshold. Hence if no bound state of H^d exists no bound state of Hint exists. It is usual to 
adjust the potential in an essentially one particle Hamiltonian such as H^d so that it tends 
to zero as ri tends to infinity. Thus the problem of showing that H^d and hence ffint have 
no bound states reduces to the problem of showing that the potential, V{r), cannot support 
a bound state, where 

V{r) = ^ + Eo{r)-ES\ (3.16) 

and 

E^^^ = lim Eo{r) = ->2 , (3.17) 

r— >oo 

is the ground-state energy of H. 

V{r) is a central potential, i.e., it is spherically symmetric. Though it depends on the 
single radial variable, r, it is a potential in three dimensions. This is very important. Any 
attractive we ll p otential can bind a particle in one dimension (see, for example. Landau 
and Lifshitz [22|) but in three dimensions it has to exceed a critical strength to produce 
binding (see, for example, Dyson on p. 1225 of 0, Wu and Ohmura HH, and Sec. fFTnl of 
Appendix). 

Ways of calcula ting the number of bound states which V{r) can support have been 



extensively studied [8^|83[. As V{r) more rapidly than r ^ as r ^ 00 and behaves like 
as r — > 0, a necessary but not sufficient condition for V{r) to be able to support A^^''* 
bound states corresponding to angular momentum I for a particle of mass m is that 

< ' / r[-y_(r)]dr , (3.18) 



{21 + 1)^2 
where 

VI (r) = V{r) if V{r) < 



y_(r) =0 if V{r) > . ^^'^^^ 



This is usually referred to as the Bargmann-Schwinger result [78|. Also important informa- 
tion about the number of bound states can be obtained by an analysis of the phase shift of 
particles of mass m scattered by V{r) (see, for example, Ref. S^l). 



The exact number of bound states of the system for a given / value can be obtained by 
determining the number of zeros (other than at r = 0) in the regular solution to the radial 
Schrodinger equation with potential V{r) and energy E = [S^. This is usually done by 
step-by-step numerical integration on a computer. 

Unfortunately, Gertler et al. [85] found that when m is the mass of the positron, V{r) as 
given in Eq. ()3.16|) could support one bound state. Thus their attempt to prove that e"'"H 
has no bound states failed. 

Fortunately, Aronson et al. js^ were able to get round this difficulty. It is reasonable to 
assume that the ground state of (e"*", H) will be an s state. For such a state, in the infinite 
proton mass approximation H^^t is of the form, 

Ant = t(ri) + t(r2) + ( + ) C + ----— , (3.20) 

\2mirf 2m2r|/ ri r2 ri2 



where 



t{r,) 



2mi rf dvi \ * dr. 



sine 89 \ dO 



(3.21) 



and 9 is the angle between the directions of the vectors ri and r2- 

In this case if we make the adiabatic approximation of fixing ri, not ri as previously, we 
obtain the adiabatic Hamiltonian 



p ' 



(3.22) 



where 



H' = tin) + - + E',{n) 



n 



(3.23) 



/(r2, 9) is a unit operator similar to /(t'2) and Eq^vi) is the ground-state eigenvalue of 



1 



1 



2ra\ri 2ra2r2 



r2 ri2 



(3.24) 



As 



(^|£^|^)>0, (3.25) 
for any allowed square- integrable function, \l'(r"i, r2), it follows that for functions of s sym- 



metry 



However, it can also be shown 



that for functions of S symmetry 



(3.26) 



(3.27) 



In addition, H'^^ and ^^int have the same continuum threshold, 
bound states, this is also true for Hint- 

The potential V'{r) associated with H'^^ is of the form 



im2. Thus if i^'j has no 



V'{r) = - + E',{r)-E'^ 



(0) 



(3.28) 



Aronson et al. showed that it is extremely unlikely that V'{r) can support a bound state 
and hence it is highly probable that no bound state of (e"*", H) exists. This result is supported 
by information from scattering calculations [84. ..87l| . Aronson et al. showed further that it 
is highly probable that no bound state of e"'"H exists for nii < 1.46 me and m2 = rrie. 

However, Aronson et al. 8^ were unable to establish their conclusion rigorously. Ar- 
mour js^ made their method of proof rigorous. He first of all calculated a very accurate 
wavefunction for the system described by H'^ using the variational method and basis func- 
tions in terms of prolate spheroidal coordinates. He was able to calculate a good lower 
bound to V^'(V) using this wavefunction and the method of Temple and Kato (see, for in- 
stance, Ref. jSSi) and show that this lowest bound, and hence V'{r), could not support a 
bound state. In a later paper. Armour and Schrader |90| showed that if mi = 1.51 me and 
m2 = me, no bound state of the e"'"H(e"'", e~, poo) system exists. To date this is the best 
lower bound on the critical positron mass required for binding. 



So far we have assumed that the mass, m^, of the proton in the (e"*^, H) system is infinite. 
Suppose we choose units in this case so that particle 2 has unit mass. Then 

Hint = -T^'^l, - ^V^, + ---- — , (3.29) 

2mi ^ 2 ^ ri r2 ri2 

where rhi = mi/m2. We can see from this that it is the ratio mi/rm2 which determines 
whether or not a bound state exists. It follows that Armour's result for the case 
mi = 1712 = shows that this Hamiltonian has no bound state for mi = m2, whatever the 
value of 1712- 

It follows from this that if the system is to have a bound state when is finite and thus 

^int = -T^'^l, - T^'^l, - —v., ■ v., + ----— , (3.30) 
2/ii 2/i2 7713 ri r2 

this must be due to the presence of the mass-polarisation term. As pointed out earlier, 
(e~, e"*", e~), and hence (e+, e~, e"*"), is known to have a bound state [Ho,!!^. Its existence 
must be due to the large mass-polarisation term in this case. 

This term can be taken into account by an interesting method due to Armour j9l|. We 
are free to choose the origin of ffjnt to be at any of the three particles which make up the 
system. Thus if we choose the origin to be at the electron rather than the proton. 

Hint = -T^K - T^^l - — Vs, • V.3 - - - - + — , (3.31) 

2vi 2/^3 m2 Si S3 Si3 

where 

mim2 17131712 

S3 = -r2 , si = ri - r2 , f i = ■ , /^s = : = 1^2 ■ (3.32) 

mi + m2 m3 + m2 

Note that fi^ = fi2 on account of the symmetry of the reduced masses of particles 2 and 3 
under the operation of interchange of these particles. 

Now we can write this form of Hi^t as an operator. A, in terms of Vi and r2 i.e., with 
origin at particle 3, 

A = --^V^^ - -^V^^ - —v., • v., - - - - + — . (3.33) 
2vi 2/i2 m2 ' ri r2 ri2 

Clearly, A has no physical significance. However, Armour was able to make use of it in 
a novel way to take into account the mass-polarisation term within the framework of the 
infinite proton mass approximation. 

Suppose Hi-^i does have a bound state, <l'(ri,r2), i.e., 

^int$(ri,r2) = ^$(ri,r2) , (3.34) 

where 

E < -i/i2 = ^thr . (3.35) 
We shall assume that E is the lowest eigenvalue of i^int and $ is normalised so that 



($1$) = 1 



(3.36) 



Because A corresponds to Hint when the origin is taken to be at particle 2 rather than 
particle 3, it follows that A must have the same eigenvalue spectrum as -ffint- In particular, 
A will also have lowest eigenvalue E. It follows that 



A-Hr 



$) > . (3.37) 



It can be shown to follow from Eq. ()3.37j) by straightforward manipulation (Armour 1982) 
that ^ I I 

' < E < -ifi2 , (3.38) 



Hi 



mt 



where 



and 



Kt = -T^^r, - T^^r, + - — , 3.39 

2mi 2fi2 ' n rs ri2 

Q = ^^^i±^ . (3.40) 

m3 — 1712 

The crucial point is that -f^int represents the internal motion, in the infinite proton mass 
approximation, of a system made up of a "positron" of the usual mass but with charge Q 
and an electron of mass /i2 and the usual charge. 

Now — |/i2 is also the continuum threshold for H'^^^ and {^\Hl^^\^) is an upper bound to 
the lowest eigenvalue of H^^^. Thus it follows from the initial assumption and Eq. ()3.38j) that 
a necessary condition for Hi^t to have a bound state is that H^^^^ has a bound state. Thus if 
it can be shown that H^^^ has no bound state, it follows that Hi^t has no bound state. 

As H[^^ does not involve the mass-polarisation term, the method described earlier can be 
applied to it. In the case of (e"*", H), 

Q = 1.0011 and /is = 0.9995 me . (3.41) 

Thus it is not surprising that Armour [91'] was able to show that no bound state of (e"*", H) 
exists, even if the finite mass of the proton is taken into account. 

(/i"*", e~, e"*") is a very similar system to (p,e~,e"'"). However, the mass of the muon is 
206.8 me, which is approximately a tenth the mass of the proton. Thus the effect of the 
mass-polarisation term will be greater for this system. Armour [9^ applied his method to 
this system for which 

g = 1.01 and /i2 = 0.995 me . (3.42) 

He was able to show that no bound state of this system exists. 

He also applied it to (p, /i~, e"*"), which differs from (p,e~,e"*") in that the electron is 
replaced by the which is much more massive. The {p, fi~) atom is much more compact 
than the H atom and has a very much smaller dipole polarisability. The reason for this is 
that the polarisability this scales as m~^, where mr is the reduced mass of the atom. See, for 
example, Bransden and Joachain As pointed out earlier, this polarisability determines 
the magnitude of the attractive potential between the positron and the atom containing the 
other two particles in the asymptotic region. Thus in this region, this potential is very much 
smaller for (p, /i~,e+) than for (p, e~,e+). 

For (p,/i-,e+), 

Q = 1.254 and fi2 = 186 me . (3.43) 
Armour 92j was also able to show that no bound state of this system exists. 



F. (p,p,e ) 



(p, p, e~) is a system containing an antiproton, a proton and an electron. As we are only 
considering Coulombic forces, as pointed out in Sec. IIII A| the Hamiltonian for the systems 
we are considering is invariant under the operation of charge conjugation, i.e., change of sign 
of the charge of all particles in the system, this system is equivalent to (p, p, e"^). 

In the absence of the electron, the proton and the antiproton can form a strongly bound 
and very compact protonium atom with ground-state energy — |yUpp, where /ipp is the reduced 
mass of the protonium atom. The expectation value of the internuclear separation in this 
state is 

1^ = 1.6 X lO-^ao , (3.44) 

where ao is the Bohr radius. 

The ground-state energy of protonium is the energy of the lowest continuum threshold 
of (p, p,e~). It can be seen that in this state the expectation value of the dipole moment 
resulting from the equal and opposite unit charges is very small indeed. Furthermore, the 
probability that the dipole moment has a value greater than 0.01 a.u. is less than 10~^. 
A consequence of this is that the dipole polarisability of the (p, p) atom, that determines 
the magnitude of the attractive potential between the electron and ground-state protonium 
asymptotically, is very small indeed, much smaller, even, than in the case of (p,yU~). 

It is clear, therefore, that no bound state of (p, p, e~) exists. As we shall see in Sec. IIV Al 
the position of (p,p,e~) on the triangular plot of the stability properties of three-body 
systems is so far into the unstable region as to be almost at one of the vertices associated 
with this region. This is to be expected in view of the very small dipole moment and 
polarisability of the (p, p) atom in its ground state. 

It is of interest to consider what would happen if the nuclei were fixed. Such a system 
has been extensively studied, in particular by Fermi and Tell er 193^ . Wightman ^9^ 
et al. jo^, Mittleman and Myerscough jo^. Tur ner a nd Fox |97j . Lev y-Leblond 19^ 
Brown and Roberts [o^, Coulson and Walmsley |lOOj . Crawford 101 1, and Turner 

It is well known that no bound state of this system exists if the dipole moment of the sys- 
tem is less than 0.639 a.u. As the nuclei have unit charge, this corresponds to an internuclear 
distance, R, of 0.639 ao. 

The Schrodinger equation fo r th is system is separable in prolate spheroidal coordi- 
nates [1^. It has been shown |lOl| that if i? > 0.639 ao, the system has an infinity of 
bound states. As is to be expected, the infinity is countable as it can be set in a one to one 
relationship with the natural numbers. In these states the electron or positron has zero an- 
gular momentum about the internuclear axis and no nodes in the part of the wave function 
that is dependent on the hyperbolic coordinate. As R is increased, further infinities of bound 
states appear corresponding to higher values of one or both of these quantum numbers, the 
first appearing when R > 3.792 ao. 




G. Muonic molecular ions 

Muonic molecular ions are systems such as (p, yU~,p), (p,yU~,d), (d, /i",d), etc. These 
systems are analogous to molecular ions such as H^, but they contain a negatively charged 
muon rather than an electron. There has been considerable interest in the properties of 
these systems as (d, /i~,d) and, in particular, (d, /x~,t), play a crucial role in the process 



know n as muon ca talys ed fusion. See, for example, Bhatia and Dracliman |103| . Ponomarev 



lOJ and Froelich |105|. 



In a muonic molecular ion, (a, where a and b are protons, deuterons or tritons, 
the dimensions of the molecular ion scale as me/m^ relative to (a,e~,6). Thus the nuclei 
a and h are brought very much closer together than in (a,e~,fe). Provided that {a,fi~,b) is 
in a state of zero angular momentum and a and b are not both protons, fusion occurs at a 
rate between 10^ and 10^^ s~^, depending on the nuclei involved. The proton-proton fusion 
process is very much slower as it involves the weak interaction. 

A muon has mean lifetime of 2.2 x 10~^ s before undergoing beta-decay. See, for example, 
Semat and Albright |l06j . Clearly, if a muon is to catalyse a large number of fusions, it is 
important that the whole process takes place in a much shorter time than the lifetime o f the 
muon. There seemed little pros pect of this until extensive theoretical work by Vesman jlOTj l 
and Gershtein and Ponomarev \lQ^ showed that {d,^~,d) and (d, yU^,t) can be formed by 
a resonant process. 

This resonant process can take place because of a remarkable coincidence. Both (d, fi~ , d) 
and (d, /i~,t) have weakly bound {J,v) = (1,1) excited states, i.e. states with angular 
momentum J = 1 and in the first excited vibrational state, v = 1. In the case of (d, /i~,d), 
this state has binding energy 1.95 eV, whereas in the case of (d, /i~,t) it is 0.66 eV. These 
binding energies are sufficiently small for energy conservation requirements to be satisfied 
by a process in which a low-energy (d, /x^) or (t,fi^) atom, in its ground-state, collides with 
one of the nuclei of a (DA) molecule {A = D, H or T), in its electronic and vibrational 
ground state, and attaches itself to the deuteron to form a muonic molecular complex, 
[(d, n~ ,d), a, e~, e~] or [(d, fi~ ,t), a, e~, e~], with the (d, d) or (d, t) in the (1,1) state. 

The complex containing (d, /i~,t) rapidly undergoes Auger decay to a state of (d, /i~,t) 
which has zero angular momentum. Fusion then occurs at a rate of 10^^ - 10^^ s~^. This is 
the fastest fusion rate of any muonic molecular ion. So far, fusion rates of about 150 fusions 
per muon have been observed in this process. The breakeven value for energy production is 
about 480 fusions per muon. 

The key role played by (d, /i~,t) in this process has focussed attention on the properties 
of muonic molecular ions. It is known that ions that contain a proton have two bound states, 
whereas (d, /i~,d) and (d, /i~j_t) have five bound states and (t,/i~,t) has six bound states. 
For details, see Ponomarev |l04| . 

As the mass of the muon is = 206.8 mg, the Born-Oppenheimer approximation does 
not give accurate results for these systems. Accurate results can be obtain ed starting from 
this approximation, using the adiabatic representation method |l04l llOOj . However, even 
more accurate results have been obtained using the Rayleigh- Ritz variational method, which 
was first applied to muonic molecular ions by Kolos et al. 110l |. This method has been 



extensively apphed to the calcula tion of the energies of the b ound states of (d, fi , t) b y, fo r 
example, Bhatia and Drachman 'llll], Korobov et al. fll^, Frolov l44J. H ara et al. Ill3, 
Hu jll4| . P etele nz and Smith 115], Szalewicz et al. (116i] . Kamir aura |ll7| - Alexander and 
Monkhorst |ll8l | and Zhen jllflt ]. For further details, see Froelich jlC 



IV. A GENERAL APPROACH TO THREE UNIT-CHARGE SYSTEMS 



In this section, we study the stabihty of three-unit charge systems as a function of the 
masses of the constituent particles. This problem has been addressed by mariy a uthors, in 



textbooks ^ 



286] , in re view papers , or in articles, see for instance 



12fiL 127l |. and references therein. 



tiy autno] 



121 



In the literature, results have been obtained on h ow the binding energy evolves one or two 
masses are varied, one of the precursors being jl28| . There are also attempts to parametrise 
the energy as a function of the constituent masses and use this parametrisation for a guess 
at the stability border, see for instance 5^ 129l |. Some papers contain thorough numerical 
investigations, where one or two mass ratios are varied. It is noted that this is a slightly 
different art to estimate accurately the energy of a well-bound system and to determine at 
which value of a mass ratio binding disappears. For instance, in Ref. jl2fl| . the stability is 
studied as a function of the mass ratios mi/m2 and m^/mi and an interesting comparison is 
made with existing results, but the matching between "molecular" states similar to and 
atomic states similar to H~ gives an unphysical spike in the drawing of the stability border. 
In Ref. ITii . an astute changes of variables makes it possible to use harmonic-oscillator 
type of wave functions, but while the binding energy of symmetric states with m2 = are 
accurately computed, the stability domain of states with m2 ^ extends too far, with, 
e.g., (M"'",M~,m^) leaving stability for M/m > 2.45, while Mitroy P] found it unstable. 
Clearly, more cross-checks of the published results are needed. 



A. Triangular plot 

Consideration of the results already described shows that stability is a property of systems 
with either mi and m2 both large or close to each other. In other words, stability requires 
|m^^ — mj^l to be small compared to m]"^. This suggests that it would be advantageous to 
consider stability as a function of the inverse masses instead of the masses rrii themselves. 
This is confirmed by observing that the inverse masses enter the Hamiltonian linearly and 
consequently the binding energies have simpler monotonic and convexity properties in terms 
of these variables. 

One can combine inverse masses and scaling and represent the stability domain with the 
normalised coordinates 

Q!i = ^ h ZT ) ai + a2 + a3 = l. (4.1) 

rrii + m2 + 




FIG. 1: The domain of possible inverse masses ai, such that ^ a, = 1, is an equilateral triangle. 

As seen in Fig. ^ each system can be represented as a point inside an equilateral triangle 
A1A2A3, the inverse mass ai being the distance to the side opposite to Aj. This is equivalent 
to barycentric coordinates. 

In this representation, the shape of the stability domain is shown in Fig. El 



Ai 




FIG. 2: Schematic shape of the stability domain in the triangle of normalised inverse masses for 
three unit-charges (+1, —1, —1). 



B. General properties of the stability domain 

The patterns of the stability and instability regions shown in Fig.|21result from three main 
properties, besides the obvious left-right symmetry corresponding to 02 as exchange. 
They are: 

1. The domain of stability includes the symmetry axis where a2 = 03. 
To our knowledge this result was first pointed out by Hill [20], who used the variational 
wave-function 

= exp(-ari2 - 6^3) + {a ^ b) , (4.2) 

already used by Hylleraas j^H for demonstrating the stability of H~ and Ps~. If the 
Rayleigh-Ritz principle is combined with the virial theorem (see Sec. IA5|) . then only the 



ratio h/a has to be adiust ed. T he apphcability of the virial theorem to variational solutions 
was noticed by Hylleraas |l3Clj and Fock |l3l| . 

2. The instability domain including A3 is star-shaped with respect to A3. 
The same holds of course for A2. We have seen that A3 does not correspond to a stable 
configuration, for a point-like protonium atom does not bind an electron. Imagine a straight 
line from A3 toward the inner part of the triangle, as pictured in Fig. IHl Moving along this 
line means keeping the mass ratio ai/a2 constant while decreases. A suitable rescaling 
provides us with a system of inverse masses 

^ (4.3) 



«1 + Ct2 ' «! + ^2 ' 1 — «3 



In this rescaled system, the masses mi and m2 are both constant, and thus the threshold 
energy E^'^^ is fixed, while increases, thus strengthening the binding of any 3-body bound 
state. This means that once the stability domain is entered, the system becomes more and 
more strongly bound, as long as the (1, 2) atom is the lowest threshold. 




FIG. 3: Star-shape behaviour: the line A3M shown crosses the stabihty border only once. This is 
because the binding energy with respect to the threshold increases monotonically on progression 
from the border to M. 

3. Each instability domain is convex. 
Imagine two points M' = {a'j} and M" = {a'-} both corresponding to instability with 
respect to the same lowest threshold, say the (1,2) atom. The same rescaling as above 

~ ry ■ 

M = {a,} M = {A}, A = (4.4) 

«! + a2 

transforms the half triangle where ^ a, = 1 and 02 < c^s into a triangle in which f3i+ (32 = 1 
is fixed. If one moves on the straight line segment between M' and M", the images of M' 
and M", say 

M = xM' + (1 - x)M" , (4.5) 

one deals with an Hamiltonian 

H{x) = H{0) + X [m - P'DpI + {(3', - (3'^)pI] , (4.6) 

where x enters linearly. Thus its ground state E{x) is a concave function of x (see Appendix). 
If E{0) > and E{1) > E^^\ then E{x) > E^^) if < x < 1. This means this instability 



domain is convex in the plane. Now the convex projection M M which moves the 
system back to the plane transforms any convex domain into another convex domain. 
The convexity property is illustrated in Fig. 0] 




FIG. 4: Convexity: if the points M and N both lie on the same stability border, any intermediate 
point on the straight segment MN belongs to an instability domain. 



C. Shape and width of the stability band 

is markedly more stable than H~. A measure of this is the dimensionless ratio 

9 = , (4.7) 

that compares the three-body energy Eq to that of the lowest threshold £'th- Values of g 
are available in the literature (see the many papers cited in the previous and the present 
sections) . They are set out in Table HI for symmetric configurations. 



TABLE I: Relative excess energy for some symmetric configurations (M^,m^,m^) 

State M/m g 

H+ 1836.15 0.19495 

(;u,d,d) 17.666 0.122 

(^,p,p) 8.8802 0.100 

Ps" 1 0.047982 

(^,e,e) 0.483793 0.05519 

H- 0.0005446 0.0553 

ooH- 0.0555 



The corresponding curve is plotted in Fig. El 




FIG. 5: The function g{a) shows the relative excess of binding, relative to the threshold, for 
{M'^ ,m~ ,m—), as a function of a = / [M^^ + 2m^^). 

The numerical calculations used to calculate g should be very accurate as they have been 
checked by several authors. It is interesting that g{a) is not minimal for a = 0, which 
would correspond to ooH~. The minimum can be estimated to be located at about a = 0.23, 
intermediate between H~ {a ~ 0) and Ps~ {a = 1/3). This explains why the stability band 
is narrower in this region. 

The link between the relative excess of binding, g, and the width of the stability band, 
6, can be made more precise. The decomposition described in Appendix, (Sec. lA 8|) reads 

H{ai, a2, as) = H{ai, a23, a23) + HpI - Ps) > (4-8) 

with = (ct2 + a3)/2 = (1 — ai)/2 and A = (02 — a3)/4. As a function of A, for given 
«! (and hence 023), the ground-state energy E is concave, and maximum at A = 0. It thus 
fulfils 

E{ai, a2, as) < E{ai, 023, 023) = [1 + fl'(«i)] -£^th(ai, "23, ^23) • (4.9) 
But the threshold energies are exactly known. Thus the above inequality reads 

1 + tti — (0:3 — 0^2) 

E(ai,a2,tt3) < Eti,{ai,a2,a3)[l + g{ai)] — . (4.10) 

1 + ai 



Hence stability is guaranteed if the last factor is less than 1, i.e. 




(4.11) 



FIG. 6: Definition of the width of the stability band. 



D. Applications of the convexity properties 



To illustrate the potential of the stability plot, let us take the example of the (1^, 1^) 
configurations. For m = 1, we have the positronium ion (or its conjugate), which is known 
to be stable. For m large, we have (e+,p,e~) which is unstable, as reviewed in Sec. IIIII It 
was stated there that the best rigorous limit for instability for this system is m < 1.51. 

Consider now an indirect approach. ooH~ is weakly bound. If one changes the mass of 
one of the electrons by about 10%, stability is lost. This corresponds to the point a in 
Fig. [7| A very conservative but rigorous estimate by Glaser et al. [13^ is that the system is 
unstable for mi = oo, = 1 and m2 > 1.57. 

Also, the system (1+, oo~, M~) is stable for M ^ oo (conjugate of ooH^), and unstable 
for M — i> 1 (this is again (e+, poo, e~)). The point /? in Fig. [7| corresponds to the critical 
mass M at which instability occurs. It has been studied by several authors. See Sec. IIIII 
For M = 1.51 or larger, the system is unstable. This value yields a conservative estimate 
for (3 in Fig. [3 Thus 7 at the intersection of ajS with the A2H axis (mi = 1713) is in a region 
of instability. It follows from this that (e~, m"*", e"*") is unstable at least for m > 4.6 me ^. 

This is much better than the results obtained by direct study of (e~, m"*", e"*"). 



Al 




FIG. 7: A safe limit on a and /3 ensure that 7, at the intersection of a(3 and A2H (mi = ma), and 
corresponding to a configuration (l+,m+, 1~), belongs to an instability domain. 

A similar reasoning shows that 7' in Fig. [3 belongs to the other instability region. Hence 
(e~, m"*", e"*") is unstable for at least m < 0.44 me. 

These values compare well with the numerical study by Mitroy j^, who estimated that 
(e~, m"*", e"*") is stable in the region 

0.69778 < m/me < 1.6343 . (4.12) 

If one takes for a and (3 values corresponding to m = 1.1 me and M = 2.20 mg, then one 
gets for 7 and 7' estimates corresponding to 

m/me = 0.64 (7') , 2.02 (7) . (4.13) 

The curvature of the stability border can be determined by comparing the values in ()4.13|1 
with those in ^J^. 



^ This corrects a misprint in Ref. [l2. 



E. Stability domain for excited states 



One needs to distinguish between two types of excited states. 

1. States with the same quantum numbers as the ground state 

2. First state with a (conserved) quantum number different from the value in the ground 

state 

As explained in Sec. m the relevant threshold should be identified. It might consist of an 
excited atom for states of the second category, if radiative processes are neglected. 

For this category of states, the variational principle (jAljl is almost unchanged. We need 
only note that the trial wave function used must belong to the subspace of the Hilbert 
space of the system with the relevant quantum number, negative parity, for instance. Most 
general results on the stability domain are immediately applicable for the ground state with 
a specific quantum number, as, e.g., total angular momentum £ = 2, or negative parity. 
However, there is a crucial exception: stability is no longer guaranteed along the symmetry 
axis. 

This corresponds to well-known properties: there are many excited states for without 
analogue in the case of Ps~ or H^. The stability domain is then a small island near the 
upper vertex Ai, as schematically pictured in Fig. |H1 A similar shape will be encountered in 
the next section, when we consider the ground state with a reduced central charge qi. 



Ai 




FIG. 8: Schematic shape of the stability domain for the state of [mf ,1112 ,mf) in the triangle 
of normalised inverse masses. 

For excited states with the same quantum number as the ground state (generalising 
the radial excitations of 2-body systems), the variational principle no longer holds in its 
simplest form. See Appendix |X] Hence, the convex behaviour of each instability region 
does not necessarily apply in the case of excited states. However, the star-shape property 
remains, as it is a consequence of the fact that P3 > 0. 



V. THREE-ARBITRARY-CHARGE SYSTEMS 



Once the domain of masses or inverse masses insuring stability of three unit charges 
(+1, —1, —1) has been estabhshed, it is interesting to go on and let the charges vary. This is 
to some extent an academic exercise, to better analyse the binding mechanisms. However, 
non-integer effective charges are at work in media, due to screening or anti-screening effects. 
Changing the charges may also be a very useful test of the accuracy of different calculational 
methods and can help to understand the correlation effects. 

The s tability of thre e-particle systems with arbitrary masses and charges was the subject 



of Refs. ^M, Esi, Il35j . We shall first review studies carried out on specific configurations. 



A. Specific configurations 

Detailed studies are available in the literature on how H~ and H2 survive if the charge Z of 
the proton is increased. They will be summarise before other mass and charge configuration 
are considered. 



1. H -like states 

In the limit where rUp is infinite, the Hamiltonian can be rescaled as 

// = ^-l + ^-l + ^, (5.1) 

and the last term treated as a perturbation around the unperturbed (Is)^ ground state. The 
series of perturbation converges towards the bound-state energy and wav e function pro vided 
1/Z < 1.09766, according to the detailed studies by Baker et al. )l36j . and Ivanov [ji 
This means that {00^ ,m~ ,m~) is bound for Z > 0.911. 



2. R+-like states 

Another well document system is {M^ , , m~) with a charge Z assigned to the nuclei. 
In the limit M/m — > 00, binding remains up to Z ~ 1.23667. See, for instance, |il25. 13^ 
and references therein. 

Note that for larger Z, till Z ~ 1.439, the Born-Oppenheimer potential between the two 
protons has a local minimum,but the values of this minimum lies above the threshold energy. 
Thus the ion has classical stability. In the following, we shall concentrate on genuinely stable 
states. 



B. General properties 

Let {mi, gi}, {m2, —q2\, and {m^, be the masses and charges, with, say, all qi > 0. 
By scaling, the properties of the (1, 2, 3) system depend on two mass ratios and two charge 
ratios. One can for instance choose gi = 1 and let q2 and q^ vary. 



For a given set of charges {qi}, one can search the stabihty and instabihty domains 
inside a triangle of normahsed inverse masses Oj, as defined in Eq. ()4.1|) . The same general 
properties hold, namely 

• every instability domain is convex 

• every instability domain is star-shaped with respect the vertices A2 and A3, as defined 
in Fig. ^ 

Likely to be lost are the properties i) that the frontier is symmetric with respect to the 
vertical axis (if m2 7^ ^3), and ii) that there is always stability along this axis. 

In fact, if m2 7^ rris, the vertical axis plays no particular role, and one should instead 
concentrate on the line 

(T) : (gig2)'(l - ^2) = (gig3)(l - as) • (5.2) 

This is a straight line issued from A'^^, the symmetric of Ai with respect to A2A3, corre- 
sponding to —ai = 0:2 = as = 1. Some examples are drawn in Fig. IHl 




FIG. 9: The threshold separation (T) is the line where the (1,2) and (1,3) atoms have the same 
energy. It is drawn here for g2/<?3 = 1 (symmetry axis), 1.5, 2, 3 and 4. 



In the {q2, qs} plane, for gi = 1 and some given masses, resulting into reduced masses 
/ii2 and /X13, (T) is a line issued form the origin. Examples are drawn in Fig. HOI 



FIG. 10: The threshold separation (T) is the Une where the (1,2) and (1,3) atoms have same energy. 
It is drawn here for a ratio of reduced masses /U13//U12 = 1 (symmetry axis), 2, 5, and 10. 

There is always a concentration of stability along this threshold separation (T). This can 
understood as follows: near (T), one can seek a wave function of the form 

^ = {[1,2],3} + {[1,3],2}, (5.3) 

where each term represents an atom linked to the third particle. Along (T), the two term 
have comparable energy, and can interfere maximally. The triangular diagram has different 
possible shapes. Let us assume q2 < qs without loss of generality. 

• If Q'2 < Q'a < 1, all niass configurations correspond to a stable ion. 

• If 5'2 < 1 < Q'3, the entire region on the left of (T), where (1,2) is the lowest threshold 
give stable configuration. There is instability at least near A3. 

• If 1 < gi < gs, there is instability at least near A2 and near A3. 

The basic observation is that if, for instance, q2 < 1, the (1,2) ion is negatively charged. 
Particle 3 feels a potential that is asymptotically a Coulomb attraction, producing bound 
states below the (1,2) energy. This state of affairs holds in particular in the limit where 
nil = nT'2 = 00 (vertex A3), and thus stability extends to the entire region from A3 to (T), 
due to the star-shape property. 

A plausible conjecture for the stability/instability frontier is shown in Fig. El for charge 
configurations (1, —q, —q) with 2 3 symmetry. For q < 1, there is stability everywhere. 
For q = 1, we are in the case studied in Sec. lIVl For q slightly larger than 1, the stability band 
shrinks. It presumably breaks down somewhere between Ps~ and H~, where the fraction of 
binding is smallest (see discussion in Sec. IIV Then stability around H~ disappears, and 
stable ions are found only in the region. Stability there is lost where q becomes larger 
than 1.24, as seen earlier in this section. 




^ < 1 



q = l 



9>1 




9 ~ 1.09 1.1<Q<1.2 9 > 1.24 

FIG. 11: Shape of the stability and instabihty regions for charge configurations (1, — g, —q). 

Some guess at the triangle plots are proposed in Fig. El for asymmetric configurations. 




93 < 92 < 1 




93 < 92 = 1 

^3 not too small 




93 < 92 = 1 
§3 small 



93 < 1 < 92 

93/92 not too small 



93 < 1 < 92 = 1 
^3 small 



FIG. 12: Shape of the stability and instability regions for some asymmetric charge configurations 
(1,-92,-93)- 



These stability domains would require thorough numerical investigation or more rigorous 
results. For instance, it has never been established that the domaii i som etimes splits into 
two parts, as suggested by one of these figures. What can be proved 134l | is that the border 
shrinks if qi decreases (for the other gj fixed) and moves to the right as q2 increases. 



Another point of view, or, say, cross section of the parameter space consists of looking 
at stabihty for given masses {rrii}, as a function of the charges. In this plot, the threshold 
separation (T) is a straight line starting at the origin. 

Suppose, for instance, that q2 increases. In the region where (1, 2) is the lowest threshold, 
one can rescale all charges by a factor l/y/q2, so that the attraction qiq2 (and hence the 
threshold energy), and the repulsion strength ^293 remain constant. The attractive term qiq^ 
decreases, so stability deteriorates. If both 52 and q^ are increased by a factor A, a similar 
rescaling leaves the attractive terms unchanged an d inc reases the repulsion by A. Again, 
one moves away from stability. It has been shown |l35j that each domain of instability is 
convex in the {l/g'2, l/^s} plane (with, say, gi = 1 fixed.). This is illustrated in Fig. IT^ 




FIG. 13: Shape of the stability domain in the plane of inverse charges (g^ ), with normalisation 
Qi = 1. The lines = 1 and = 1 are either asymptotes, or part of the border, starting from 
a value of or ^3^^ which might be less than 1, unlike the case shown in this figure. 



C. Critical charge for binding two electrons 

The s ymm etric case, where m2 = = m, and qi = 1, q2 = q3 = Z , has been studied by 
Rebane 125l |. who addressed the following question: what is the maximal value, Z^r, of Z 



allowed for which the system is bound? This is equivalent to finding the minimal charge qi 
required to bind two identical unit charges. 

It has been seen that q ~ 1.098 for mi/m 00, and q ~ 1.24 for nii/m —>■ 0. The results 
for the general case are given in Fig. 




0.2 0.4 0.6 0. 



FIG. 14: Estimate of the maximal value Z^r of the charge Z allowed to bind i {m\ ^ 1^^, a 
function of the normalised inver se m ass a\ = 1/(1 + 2mi). The curves are rough fits to the upper 
and lower bounds given in Ref. |l25l |. 



There is a striking analogy between this curve and that of Fig. |S] showing the relative 
binding energy for Z = 1. Again, the minimum is not exactly in the H~ limit, but for 
cti ~ 0.2. 



D. Further limiting cases 

In addition to H~-like, H^-like and other symmetric cases, results have been obtained for 
the following configurations: 



1. Nearly symmetric states 

Equation ()4.8p was a decomposition into symmetric and antisymmetric parts for q2 = 
and 1712 slightly different from m^. Similarly, one can consider the case m2 = m^, i.e., the 
vertical axis of the triangle, and q2 and q^ slightly different. Introducing the strength factors 
■51 = q2(l3, S2 = qsqi, and S3 = qiq2, and denoting E[si, S2, S3] the ground state energy, one 
gets by the same reasoning 



E[si, S2, S3] < E 



S2 + S3 S2 + S3 



(5.4) 



2 ' 2 

with the later energy being [l + g) times its threshold energy, as per Eq. ()4.7p . If for instance 

Si = l, S2 = l-X, S3 = l+X, (5.5) 



stabihty is ensured if 



:i + xY<l + g 



(5.6) 



which can be solved for the charges, rescaled to gi = 1, with the conclusion that stability 
remains at least up to the configuration 




FIG. 16: Schematic shape of the stabihty domain in the Born-Oppenheimer hmit. The heavy 
particles have charges q2 and ^3. The charge of the hght particles is set to gi = — 1. 

When «! varies, the minimum of g{ai) is between 4% and 5%. The, if x < 0.02, the vertical 
axis of the triangle belongs to the stability domain, as shown in Fig. El 



Ai (T) 




FIG. 15: For charges qi = 1, 52 = 1.02, and 53 = 0.98, the stabihty domain extends at least from 
A2 to the median axis of the triangle. 



2. The asymmetric Born-Oppenheimer limit 

Consider the stability domain in the {q2, qs) plane in the limit where both m2 00 and 
m3 —* 00. The domain includes the unit s and square and extends till ^2 = 93 — 1-24 on the 
symmetry axis. This is ill ustrated in Fig. 1161 

It can be shown 134l | that the boundary leaves the square vertically at the point (^2 = 
1) 93 = 0) with a behaviour 

* - 1 " ''j^^ ■ 

A similar pattern is observed near (g2 = 0, gs = 1). 



E. Numerical investigations 



Modern computers and new computational algorithms make it possible to perform intri- 
cate few-body calculations in a very short time, and thus to repeat the calculations with 
changes in the paramet ers. 

For instance, Krikeb |l39| has drawn the stability area (more precisely its minirnal exten- 
sion, as his calculation was variational) for different choices of constituent masses |l35L Il39j | . 
Some of his plots are reproduced in Figs. 117^001 There is always a pronounced spike near 
the threshold separation (T). 




FIG. 17: Variational estimate of the stability domain for mi = oo and ■m2/m^ = 1, full view (left) 
and enlargement of the spike (right). 
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FIG. 19: Same as Fig. \17\ for m2/m^ = 1.5. 
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FIG. 20: Same as Fig. HH for 7712/^3 = 2. 



F. Open issues 

Some questions remain, for instance: 

• In the (g2, Q'3) plot, are there always points for which ^3 > 1? 

• What is the overall shape of the areas of stability? This means a curve /(^2,Q'3) = 
such that beyond this curve, stability never occurs, whatever values are assumed for 
the constituent masses. A guess is proposed in Fig. EU 




FIG. 21: Variational estimate of the domain in the {q2,Q3) plane, for qi = 1, for which stability 
occurs at least for some specific values of the mass ratios. 



VI. FOUR-CHARGE SYSTEMS 



In this section, we consider systems [mf, m^, , m^) with arbitrary masses rrii and unit 
charges gi = g2 = +1 and qs = q4, = —I- The question we address is: for which values of 
the masses rrii is the system stable against spontaneous dissociation? The answer, of course, 
only depends on the mass ratios and not on the masses themselves. 



A. Lowest threshold 



The first stage is to identify the lowest threshold. Without loss of generality, one can 
assume mi > m2 and > m^, so that the lowest threshold consists of two neutral atoms 
is (m^,mj) + (m^,mj). This can be seen directly from t he B ohr formula. The resu lts is 
in fact more general, as noticed by Bertlmann and Martin 14Cl| | and by Nussinov |l4l| : if a 



given two-body potential provides two constituent of masses m and m' with a ground-state 
energy E2{m,m'), then 

£2(7711,7713) + E2{7n2, 7714) < E2{7ni,77l4) + E2{77l2, 771^) , (6.1) 

when mi > m2 and ms > 7714. This is easily shown using the property that the two- 
body energy i?2(m, m') is a convex function of the inverse reduced mass m~^ + m'~^.. For 
instance, in the limit of a strictly flavour-independent confining force between a quark and 
an antiquark, one anticipates a mass inequality 

(QQ) + (qq) < 2(Qq) , (6.2) 



which is observed in the actual spectrum of mesons |14^ 

There is another type of threshold with a three-body ion and an isolated charge. For 
instance, (m^m^mj) + mj is certainly the lowest threshold if mi = m2 = ma ^ m4, as 
typical Coulomb energies grow proportionally with masses. However, an ion with charge 
±1 can always bind a charge =f1, as their interaction reduces to a Coulombic attraction at 
large separation (see Sec. lA iJ of Appendix El) • Thus, if the lowest threshold is an ion and 
a charge, the system is stable. 

In short, to establish stability, it is sufficient to get below the lowest atom-atom threshold, 
i.e., to find a trial wave function with a variational energy -E4[\E'] such that 

E^l^] < E2i7ni, mg) + ^2(^2, 771^) . (6.3) 



Nevertheless, it is interesting to determine which is the lowest threshold. In Ref. |14c 
for instance, dealing with hydrogen-like configurations (M+, M+, m~, m~), it is shown rigor- 
ously that the lowest threshold always consists of two atoms, namely, in an obvious notation, 

2E2(M,m) <E3(M,M,m) , 

2E2{M,m) < E3{m,m,M) . ^ ' ' 

In practical calculations, knowing that the lowest threshold of a system is made of cluster 
Ci and cluster C2 suggests a possible trial wave function associating the internal motion 
within Ci and C2 with the motion of Ci relative to C2. 



B. Specific systems 



1. H2 

The binding mechanism of the hydrogen molecule, H2, is described in several textbooks. 
A b eauti ful picture is provided by the method of Born, Oppenheimer, Heitler and Lon- 
don |l44| . In the limit where the proton mass M is very large compared with the electron 
mass m, the protons move almost classically in an effective potential V{R). Besides the direct 
electrostatic repulsion —l/R, this potential is generated by the energy of the two electrons 
in the static field of the two protons separated by a distance R. At very large R, the system 
consists of two independent atoms, and one fixes the potential so that V{oo) = 0. At very 
small R, the repulsion dominates, and V +00. In between, there is a well-pronounced 
pocket of attraction with V < 0. In the extreme limit where M ^ 00, the minimum of V is 
the binding energy of the molecule. This minimum is deep enough to guarantee that for the 
actual molecule, binding survives numerical uncertainties, finite proton mass and quantum 
effects. 

In fact, the potential V between the two protons is deep enough not only to bind the 



molecule, but also to allows for many excited states. See, e.g., |145 . 

Several accurate, r ion-a diabatic, calculations of the H2 energy and properties have been 
pubhshed, see, e.g.. Note that it was once believed that a problem existed, with 



a discrepancy between the experimental ionisation energy and the theoretical one |14 
Nowadays, the problem has disappeared. 



2. PS2 

The stability of (e+, e+, e~, e~) with respect to dissociation into two positronium ions 
(e"*", e~}_was speculated by Wheeler A first study by Ore in 1946 led him to con- 

clude |l48l | "with reasonable assurance that this structure is not stable against disintegra- 
tion into two bi-electrons." He used, however, a sharply truncated harmonic-oscillator basis 
which is not appropriate for a weakly-bo und molecule. 

The following year, Hylleraas and Ore |l49l | introduced a more suitable trial wave-function 
and the stability of the positronium molecule {m~^ ,m~^ ,m~ ,m~) was established, and con- 
firmed in several further studies using more refined wave functions. In natural units, the 
positronium atom (m"'",m~) has a binding E2 = —1/4, so that the dissociation threshold 
for the molecule is Ett = 2E2 = — 1/2. It is convenient to define the fraction x of binding 
below threshold as 

E(e+,e+,e-,e-) = (l + x)Eth . (6.5) 



The result of Ref. |149| . corresponding to a; ~ 0.0 084, has been improved b y sev eral authors. 
The lates t (and best) values are x = 0.03186 and x = 0.03201 ^E^. See, also, 

iHlili. 



The proof of stability of Ps2 by Hylleraas and Ore |149| relies on an elegant variational 
method. They first got rid of the scale by noticing that if \E'(rj) is a trial wave function, 
with norm and expectation values of the kinetic and potential energies written as 

n=(^|^), t = (^|T|^), v={^\V\^), (6.6) 



then the best rescahng of the type = \E'(rj/A) yields a minimum 



E = -^, (6.7) 

which corresponds to {0|T|(/)) = —{(j)\V\(j))/2, i.e., the same sharing of the kinetic and po- 
tential energies as for the exact solution. This extension of the virial theorem to variational 
approximations is well known jlSOl Il3ll | . 

The frozen-scale wave function of Hylleraas and Ore only contains a single parameter: 

= exp[-(ri3 + ri4 + + ^^24/2)] cosh[/?(ri3 - - + ?^24)/2] . (6.8) 

Explicit integration leads to (a misprint in Ref. |l49l | is corrected below) 

_ 33 33 - 22/32 + 5/3^ _ 21 _ 3/3^ 21 - 6/3^ + (3^ 

~ 16 ^ 16(1-/32)3 ' 8(1-/32)3 ' ^ ' 



19 21 - I8/32 + 5/3^ 1 



6 4(1 -/32)3 (1 -/92)2 



5/32 1 7 (l-/32)\ 



4/?4 8/?2 4/?6 1-/32 



to be inserted in ()6.6|) and ()6.7|) . leading to a minimum ii^ = —0.5042 near (3"^ = 0.48. 
More details of the derivation of Eg. (16.911 are provided in Appendix El It is rather easy to 
generalise the calculation of Ref. |l49| to a trial wave function of the type 

^ = ^Ciexp[-ai(ri3 + ru + ^23 + ^24)] cosh[6i(ri3 - ru - r-23 + ^24)] , (6.10) 



but one does not gain much 154l |. As analysed for instance by Ho j69j and by the authors he 



quotes, some explicit ri2 and r^^ dependence is needed in the wave function to ir nprove th e 
accuracy. Such a dependence is in clude d in all accurate variational calculations |l50L Il5lj , 
as well as in Monte-Carlo variants |l52| . 

As Ps2 is weakly bound, one might guess that there are no bound excited states, and, 
indeed, for many years, no excit ed state was found or even searched for. However, a recent 



highly accurate calculation jl51| indicates the existence of a state Psj with angular momen- 
tum L = 1 and negative parity. This cannot decay into two positronium atoms in their 
ground state, as one of the atoms should be in a state with L = 1. The threshold energy for 
spontaneous dissociation is thus 

^t*h = (1 + ^) = -0.3125 , (6.11) 

while Psg is found |l5l| to have an energy E* = —0.3344, i.e., bound by a fraction x* = 0.07, 
where x* is defined as in ()6.5|1 . 

Electromagnetic transitions from Psg to Ps2 might be used as a signature of the formation 
of positronium molecules. The main decay channels, Ps2 — > 47 or Psg — > 47, are not of great 
use, unfortunately, as the energy of the emitted 7-ray differs very little from those coming 
from annihilation of ordinary Ps atoms. This is discussed in Ref. |l55j |. 



3. PsH 



The hydrogen "hydride" PsH is the system (p, e+, e^, e^), or in our other notation, 

{M^ ,m'^ ,m~ ,m~) in the hmit of a large M/m mass ratio The chemistry of positrons 

is beyond the scope of this review. See, for example, Ref. |l5fiL p. 263]. What is of inter- 
est to us is the fact that, if annihilation is neglected, there are many examples of stable 
configurations where a positron is attached to an atom or an ion . 

The stability of PsH is not completely obvious, if one starts building this configuration 
by adding the constituents one by one. In a possible such construction scheme, one is faced 
with the instability of (p, e+, e~), as seen in Sec. IIII E| it is thus a pleasant surprise that the 
whole system is stabilised by adding a second electron. On the other hand, when one starts 
from a Ps~ ion and a proton, one expects some attraction between them. The question is 
whether this attraction brings the system down not only below the (Ps~,p) threshold but 
also the (H, Ps) threshold. A similar question arises when the system is viewed as a stable 
H~ ion combined to a positron. 

The first proof of stability of PsH is due to Ore jl57l |. He used a trial wave function 
similar to ()6.8|) . namely, with a similar notation (1 is the proton, 2 the positron) 

^ = exp [-k{n3 + aru + Pr24)] + {3^4}. (6.12) 

Again, the scale parameter k can be fixed by the virial theorem, and the minimisation is 
restricted to the parameters a and (3. Some details on the calculation of the relevant matrix 
elements are given in Appendix O For a = 0.25 and (3 = 0.5, Ore obtained an energy 
E = —0.75256 corresponding to a fraction x ~ 0.00334 of extra binding, defined as in 
Eq. ()6.5p . With a slightly more general wave function 

^ = exp [-(ri3 + 7r23 + ar^ + Pr24)] + {3^4}, (6.13) 

one obtain s a better value x = 0.0112. With a polynomial added, Navin et al. obtained 
X ~ 0.039 |l58j . 

Positronium hydride was identified experimentally in 1992 |l59| . Me anwhile, variational 



calculations have been pushed to higher accuracy. See, for instance, |150l . I155L [1 



161 

I62L I63L 164 165l . 166^ , and references therein. For instance, Zong-Chao Yan and Ho 166 



obtained x ~ 0.052. 

The system (p, e+, e~, e~) also has a resonant state with angular momentum L = 0. It is 
located aX E = —0.602 with respect to complete dissociation, and thus is relatively narrow 
(width r = 0.0028) since it is below the threshold for dissociation into a hydrogen atom and 
an excited positronium atom (E'th = —0.5625) [16^]. There is also a metastable L = 1 state 
at E = —0.615 |16l| . which is thus below the threshold Eth = —0.5625 for dissociation into 
a hydrogen atom and a Ps atom with L = 1. The situation is similar for the metastable 
L = 2 state at ii^ = —0.566 jl6ll |. which is below the threshold Eth = —0.5267 corresponding 
to H and Ps with L = 2. 



C. Properties of the stability domain 

1. Tetrahedral representation 

The stability domain of the system, (mj*", m2 , , rn'^), made up of four particles, two 
having the same positive charge and two having the opposite negative charge, interacting 



only through Coulombic forces, can b e represe nted in terms of points in the interior, or on 
the surface of a regular tetrahedron 167 . 168j |. This is a generalisation of the equilateral 
triangle diagram in the case of three bodies, tw o of w hich have the same positive charge and 



123l |. Note that as the Coulombic potential 



the third has the opposite negative charge 
between charges is invariant under charge conjugation, stability properties are invariant 
under this operation. 

Let ai = l/rrii, were is the mass of particle i. As only Coulombic forces are being 
considered, the stability properties are dependent only on the ratios of the masses. It is thus 
possible to normalise the so that Y2t=i = 1- The sum of the distances from each face 
to any point X, inside or on the surface of the tetrahedron, is a constant. If we choose the 
length of the sides of the tetrahedron to be ^^(3/2), this constant is 1. Take aj to equal 
the distance of the point X from face j. Then the set X (ai,a2,a3,a4), with < aj < 1, 
Y2i=i Q^i = 1 represents all points within the tetrahedron and on its surface and hence all 
possible mass ratios of the system {mf, m^, mj, mj). 

The regular tetrahedral diagram is shown in Fig. |221 The vertices correspond to cases 
when three of the particles are infinitely massive. The side AB corresponds to the set of 
inverse masses (1 — x, x, 0, 0) with < a: < 1, where two of the particles, in this case 3 and 
4, are infinitely massive, and similarly for the other sides. 




X (qi, Q!2, as, cXi) 



D (0, 0, 0, 1) 
0, 0.5, 0) 



B (0, 1, 0, 0) 



C(0, 0, 1, 0) 



FIG. 22: The regular tetrahedron diagram representing the system (m^, m^, rrig , 70,4 ). Faces 
BCD, ACD, ABD and ABC are labelled 1,2,3 and 4 respectively. Points E, F, G, H, J and K are 
the midpoints of the sides of the tetrahedron. Regions Ri and R2 correspond to the domain of 
instability on the triangular surface bed. 



It is interesting to consider the following special configurations: 

• H-H corresponds to (p,e+,p,e-), i.e. to (0.4997,0.0003,0.4997,0.0003) which is very 
close to i?(0.5, 0, 0.5, 0), the midpoint of AC. Any permutation of particles 1 and 2 



and/or 3 and 4 also corresponds to H-H and thus the midpoints F,G and J of BC, BD 
and AD, respectively, also correspond very nearly to H-H. All these points are thus 
deep inside the domain of instability. 

• The plane KDC corresponds to systems with identical particles 1 and 2, with mi = m2, 
and the same charge. It can be seen that they are all stable. Similarly, the points of 
the plane ABH represent systems with identical particles 3 and 4, which are stable. 

• The midpoint H of CD lies very close to (p, p,e~,e~), which is H2 and is known to 
be bound. Similarly, the midpoint K of AB lies very close to (e"*", e"*", p, p), which 
is the antihydrogen molecule with the same bound states as H2. The segment KH 
corresponds to hydrogen-like configurations {M~^ , ,m~ ,m~). In the middle, the 
centre of the tetrahedron represents the positronium molecule. 

There are the following obvious symmetries of the stability domain of {m^, m^, , mj): 
1 ^ 2, 3 ^ 4 or {1,2} ^ {3,4} exchanges. As a consequence, the binding energy is 
stationary around the corresponding symmetry axes. For instance, it is of second order in 
Sm = 1712 — mi , f or gi ven mi -|- m2. 



Richard |167l Il68| considers that the lines AC, BC, BD and AD represent unstable 
systems whereas AB and CD represent stable systems. He gives schematic pictures of the 
stability domain on each of the faces of the tetrahedron. It takes the form of a sharp peak 
with base AB or CD, which is symmetric about the median perpendicular to AB or CD on 
the face under consideration. 

Richard describes what is known about the form of the stability domain within the tetra- 
hedron. The form on a triangular surface close to a vertex and parallel to the opposite 
face is an interesting case. This surface is in the shape of an equilateral triangle and repre- 
sents systems in which one particle is very much lighter than the other three. As Richard 
points out, in this situation the existence or otherwise of four-body bound states depends on 
whether the three very much heavier particles have a bound state that can form a charged 
'nucleus' that can bind the very much lighter particle, which is of opposite charge. It follows 
that the form of the stability domain on the above surface is the same as for three particles, 
two of which have the same positive charge and the other has the opposite negative charge 
JSS], with upper vertex on AB or CD, the sides of the tetrahedron that are considered 
to correspond to bound states. See Fig. |221 To summarise, the shape of the stability domain 
is likely to be that schematically pictured in Fig. 123 Some guessed properties remain to 
be demonstrated, and more numerical investigations would be needed for a more precise 
drawing. 
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FIG. 23: Guess at the shape of the stabihty domain for four unit charges, in the tetrahedron of 
normahsed inverse masses. 

2. Hydrogen-like configurations 

We have seen that {M~^ , ,m~ ,m~) is stable for M ^ m and for M = m, but 
this is shown with different approaches: the Heitler-London method for hydrogen and the 
Hylleraas-Ore variational method for Ps2. In fact one can show that the stability of Ps2 
implies the stability of {M^ , ,m~ ,m~) for any value of the mass ratio M/m. 

The proof follows from the property that the ground state of an even Hamiltonian in one- 
dimensional quantum mechanics is lowered by any odd admixture. If h{X) = /leven + Avodd, 
the ground-state energy, e(A), of h{X) cannot be greater than the expectation value obtained 
with the even wave-function corresponding to A = 0, which is, of course, e(0). An illustration 
is h{X) = + + \x, whose gr ound state lies at e(A) = 1 — A^/4. Here, parity is replaced 
by charge conjugation |l43l Il70| . The Hamiltonian for {M~^ , ,m~ ,m~) can be written 
as 

H =Hs + Ha , 




(6.14) 



The term is odd under charge conjugation, i.e., simultaneous Pi ^ and P2 ^ 4p4 
exchanges. It improves the binding of the system. More precisely, the ground state E 
of H is lower that the ground state of Hs alone, which is none other than the Hamilto- 
nian of a rescaled positronium molecule with an inverse constituent mass (m~^ + M~^)/2, 
which provides Hs with the same threshold as H. Thus the stability of Ps2 implies that of 
{M^ , ,m~ ,m~) for any M/m ratio. The additional binding provided by the term 
is seen on the dimensionless fraction x measuring the fraction of binding below the lowest 
threshold. In the best variational calculation, the fraction x, indeed, increases from x ~ 0.03 
for Ps2 to X ~ 0.17 for H2. 



3. Minimal extension of the stability domain 



Once the binding energy E = {1 + x)E^\^ of configurations with mi = m2 and 777,3 = ''^4 is 
taken as firmly estabhshed by detailed variational calculations, a minimal extension of the 
stability domain can be derived. 

Consider first {M~^ ,m~^ , M~ ,m~) configurations as extensions of the Ps2 case. One can 
fix the scale by setting = 1 — y and = 1 + y. Then 

HiM\m\M-,m-) = ^^^^±^±^ + V + y-{-pi + pl-pl + pl) . (6.15) 

The reasoning is the same as previously: as the last term is antisymmetric under simulta- 
neous Pi ^ P2 and P3 ^ P4 exchanges, it lowers the ground-state energy. Thus, 

Eiy)<{l + x)EttiO) . (6.16) 

Meanwhile the threshold becomes 

1 1 g.i.(0) 

^"<^) = -1(13^ - ioT^ = ■ 

Thus stability remains at least as long as < 1 — (1 + x)^^, i.e., for 

M 

0.70 < — < 1.43 , (6.18) 

m 



if on e accepts the value x = 0.032 for Ps2 |151l |. For comparison, a recent numerical study 
|l6ll | show that the actual range of stability should be very close to 

M 

0.476 < — < 2.1 , 6.19 

777 

Now, if one starts from (M"^, m+, M~, m~), and introduces four different masses rrii such 
that 

?77i > 7779 , 7777"'^ + 777^"'^ = 2M~^ , 

-1-1 -1 (6-20) 

one can rewrite the Hamiltonian as 

H{mf, m^, mj, mj) = H{M^, m^, M~, 777") 

+ \[{m^'~m^'){pl-pl) + {m,'-ml'){pl-pl)] . ^^'^^^ 

In Eq. ()6.2H) . the first term is symmetric under simultaneous Pi ^ p^ and P2 ^ P4 ex- 
changes, and the second one is antisymmetric. By the same reasoning as previously, the last 
term improves the stability of the system, without changing the threshold energy. Equa- 
tions ()6.18j) and ()6.2()|1 represent a minimal extension of the stability domain around Ps2, 
based solely on the variational pri nciple and the basic symmetries. 

One can thus look in the tables (l42| for charged hadrons or leptons which are long-lived 
and list many highly-exotic molecules which are stable. One example is the strangeness —3 
and charm +1 system {Q~ , Il~ , , p) system, as the constituent masses are (in GeV/c^) 
1.672, 1.197, 1.869 and 0.938, respectively, and satisfy Eqs. (IFTTHl) and dH^Ol)- 



Going back to more realistic cases, one can see that Eqs. ()6.18|) and ()6.20j) fail by a narrow 
margin to establish the stability of PsH from the binding energy of Ps2. 

Suppose now one assumes a binding energy E = {1 + a;(s))£'th for a hydrogen-like con- 
figuration {M~^ M~^m~m~) with mass ratio s = M/m. Introducing 



m 



rrin 



-1 



m 



-1 
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This leads to a decomposition of the Hamiltonian 

Hlm^ ,1712 ,171^ ,171^) ~ H{M~^ , ,171^ ,m^) 



5' 



[pI - pI) 



(6.22) 



(6.23) 



which implies E{mi, 1712, 1713,1714) < (1 + x{s))Et^{M, M,m,m). Since the threshold for the 
{rrii} configuration is 



Eti,{{mi}) 



-1/2) 
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we have stability at least as long as 

6 + 6' 



M-i + m- 



< 1 



M-i + in' 



1 



(6.24) 



1 + xis) 



(6.25) 



Thus if it can be shown that a system (M"*", M"*", m , m ) with mass ratio s = M/m = 2 is 
bound by a fraction x(2) larger than 1/8, then it follows that PsH is stable. 



4- Stability for equal-mass negatively-charged particles 

In the Born-Oppenheimer approximation, all molecules such as (p,d, e~,e~) 
or (p,t,e~,e~) are stable. We have seen that all hydrogen-like configurations 
{M~^ , ,m~ ,m~) are stable, even for M c:^ m. The positronium hydride (p, e"*", e~, e~) 
is stable. For M ^ m ^ M', the configuration {M^ , M'^ ,m~ ,m~) is stable, as the lowest 
threshold consists of the stable ion {M^ ,m~ ,m~) and an isolated M'~ . 

These observations make it reasonable to conjecture that all molecules (m^, m^, ,1714) 
with equal-mass electrons (ms = m^) are stable. Of course, a similar result holds for nii = 
1712- This can be understood as follows: the system (M~^ , M'~^ , in~ , in") has two degenerate 
thresholds, one with M associated with the first electron and M' with the second and another 
in which M and M' are interchanged. There are thus two ways of describing the system 
as two deformed atoms with a slow relative motion. These two components can interact 
strongly with each other, leading to binding. In the 3-body case, the same phenomenon 
leads to enhanced stability of [M^ ,m~ ,m'~) around the symmetry axis m = m' . 

The conjecture of stability for = has been checked in Ref. jl7l| . The authors 
used two variational calculations. The first one is based on the exponential wave function 
()6.13|1 . the second is the stochastic variational method with Gaussian wave functions. As all 
matrix elements are calculated in a close analytic form, there are no numerical uncertainties 
as to the values of the integrals involved, so the investigation can be considered as free from 
ambiguity. 



5. Borromean binding 



The concept of Borromean binding is inspired by the Borromean rings, which are inter- 
laced in such a subtle topological way that if any one of them is removed, the two other 
become unlocked. The name "Borromean binding" was introduced in Nuclear Physics jl72j | 
to denote those isotopes with two external neutrons which are stable against dissociation, 
while the partner with only one external neutron is unstable. The best-known example is 
the Helium family. The nucleus a = ^He is very stable and compact. The isotope ^He is 
unstable, but ^He is stable. If considered as a three-body system {a, n, n), it has the remark- 
able property that none of the two-body subsystems has a bound state. In other words, with 
short-range potentials, a three-body system can be bound, with coupling strength that are 
too weak to bind the two-body subsystems. For an inve stigation of th e domain of coupling 
constant where Borromean binding is possible, see e.g., 

Borromean binding is related to the "Thomas collapse" (l7fi| . which occurs when the rati o 
of 3-body to 2-body binding energies, E3/E2, is very large, and to the Efimov effect |l77l |. 
which is associated with the existence of very weakly bound 3-body bound states together 
with a vanishing 2-body ground-state energy E2. 

Borromean binding is defined to occur in four-body binding as the property that all 
three-body subsystems are unbound. In other words, there is no way to build the system 
by adding the constituents one by one, forming a bound state at every stage. Of course, for 
systems such as (a^ ,b'^ ,c~ ,d~), there always exist stable two-body subsystems. 

Now, Mitroy p| established that {m'~^ ,Tn~^ ,Tn~) ions remain stable for 

in' 

0.698 < — < 1.63 . (6.26) 

m 

On the other hand, Bressanini et al. |l52j found that the (M"*", m"*", M~, m~) systems are 
stable with respect to dissociation into {M^^M~) + {m^ ,m~) for 

1 M , , 

— < — < 2.2 . 6.27 

2.2 - m - ^ ' 



This result is also obtained by Varga jl78[. Thus for M/m ~ 2, the 4-body molecule is 
bound, but none of its 3-body subsystems is stable! This means, for instance, that the 
molecule (p, d, p, d) made up of a proton, a deuteron and their antiparticles, is Borromean. 
A minimal extension of the domain of Borr omean binding can be established if one knows 
the binding for mi = ms and 7712 = 777.4 |l79j . 

The reciprocal conjecture, stating that if a four-body systems has at least one stable three- 
body system, then the overall system is stable, remains, however, to be demonstrated. An 
argument is the following. Assume < and suppose that (1,2,3) is stable (if (1,2,4) 
is stable the stability of (1,2,3) is implied), i.e., with the current notation, E123 < -E13, 
and consider the system {mf , 1712 , , m^) as a function of = I/7774. For X4 = 1/777,3, 
the four-body system is stable, as seen in the previous subsection. For X4 —>■ 00, the lowest 
threshold becomes the charged (1, 2, 3) ion, which attract the charge 777^. As seen in Fig.l2H 
it is very unlikely that the curve of the four-body energy would not remain below the lowest 
threshold between X4 = l/m3 and X4 = 00. 
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FIG. 24: Guess at the shape of the four-body energy of {m\ ^m^ ,rn'^ ,m^) (dotted curve) as a 
function of X4 = l/m4, from 2:4 = l/ms to large X4. The threshold is the lowest of the (1, 3) + (2, 4) 
energy (dashed curve) and the energy of (1,2,3) (solid line), which is assumed to be stable. 



VII. SYSTEMS WITH FIVE UNIT CHARGES 



Calculating the binding energy of a five-particle system is obviously a very complicated 
task [180i . il8l[ll82j . especially when the system is loosely bound. A five-particle system can 
have various bound subsystems including 3 + 2, 4 + 1 or 2 + 2 + 1 decompositions. The 
binding energy is the energy of the system with respect to the lowest dissociation threshold. 
We shall consider the case of systems of three particles with like charges and two particles 
with the opposite charge. Systems of unit charges where four particles have identical charges 
and the last particle has an opposite charge are not expected to be bound. 

A. {m'^ ,m~^ ,m'^ ,m^) 

The simplest system consists of five particles with equal mass and unit charge. In the 
case of where the three positively-charged particles are identical fermions with spin 1/2, the 
antisymmetry requirement restricts the configuration space and no bound state exists. In 
particular, the system of three electrons and two positrons, or its conjugate, is not bound. 

The {m^ ,m~ ,m~) system is, however, bound if the three positively-charged 

particles are bosons or if one of them is distinguishable, since the Pauli principle does 
not restrict the allowed states. An example is the (e"'",e+,e~,e~,a:^) system, where 
fictitious particle which has the same mass as the electron but is distinguishable from both 
the electron and the positron. No such system exists in the real world. What is of practical 
interest is that binding often remains when some masses are changed, in particular when rrix 
is increased. For example, a system made of a Ps2 molecule and a proton, (e"'",e"'",e~,e~,p), 
is stable (see Subsecs. IVII HI and IVII II below) . 

The energies of systems of equal-mass constituents (with mass equal to rrie) are listed 
in Table |nl for both the boson and fermion cases. The Ps~ ion and Ps2 molecule are well 
known examples of such systems. 

TABLE II: Energies of A^-particle systems of unit charges and equal masses. The total charge is 
and 1 for even and odd, respectively. Atomic units are used. 



N 


fermions 


bosons 


2 


-0.250000 


-0.250000 


3 


-0.262005 


-0.262005 


4 


-0.516004 


-0.516004 


5 


unbound 


-0.556489 



In Table ITTj the energies of the bosonic and fermionic systems are equal up to A^ = 4. In 
the bosonic case, the particles are considered to be spinless and the spatial part of the wave 
function is symmetric in the coordinates of the identical particles. In the lowest-energy state 
of the fermionic system (particles with spin 1/2), the spins of the pairs of identical particles 
are coupled to zero. In this case the spin part of the wave function is antisymmetric and 
the space part therefore has to be symmetric. Consequently, both the bosonic and fermionic 
systems have symmetric spatial parts and their ground state energies are equal. 



B. {M+,M+,M+,m-,m-) 



A well-known stable Coulombic five particle system is the Hj]" molecule where the three 
protons form an equilateral triangle sharing the two electrons. The stability of this system 
is due to the slower motion of the heavier particles. If the mass of the protons were equal 
to that of the electrons, there would be no bound five particle system, as we know from the 
previous section. The question is at what heavy/light particle mass ratio is the stability 
lost in the case of {M~^,M~^,M^,m~,m~). A system of three holes of mass and two 
electrons in sem icon ductors is a realistic example of this. The stability of such systems has 
been studied in 180l |. The energy as the function of the electron/hole mass ratio is shown in 
Fig. 1231 The system is bound provided that the positively charged particles are at least five 
times heavier than the two with negative charge (0 < cr = m/M < 0.2) (m = mg, M = rrih). 
At higher value of the mass ratio, the system dissociates into a (M"*", M"*", m~ ,m~) molecule 
and an isolated ' " ' ' x , ;(M"'",m~). 

However, the ei 1. 




FIG. 25: The total energies of three-hole and two-electron systems as a function of the electron to 
hole mass ratio a = nie/mh- The energies of the relevant three- and four-body thresholds are also 
shown. 

This study shows that the Hjj' molecule would remain stable even if the protons were be 
much lighter. The system consisting of three protons and two negatively charged muons, i.e., 
(p, p, p, can be mentioned as an exotic example where a < 0.2 (see Table UTT]) . The 
mass ratio between the muon and the proton is about a = 0.11 which is much larger than 
that in the hydrogen atom {a = 0.0005). The energy of the proton-muon atom is —92.92 
a.u. The molecule formed by two proton-muon atoms is deeply bound just like the hydrogen 
molecule [any (M"*", M"*", m~, m~) system is bound irrespective of the M/m ratio, as seen in 
Sec. Ell- The binding energy divided by the reduced mass of the proton-muon atom is 0.07 
in (p, p, /i~) and 0.02 in (p, p, p, /i~). The corresponding ratios for H2 and Hj]" are 
0.16 and 0.18, that is the (p, p, p, /i~) is much more loosely bound than the H3 . The 
energy of the (p, p, /i~) ion is —102.22 a.u., corresponding to —195.14 for (p, p, /i~) -|- (p, /i~), 
to be compared to —1999.63 for the (p,p,/i~/i~) molecule, an illustration of the fact that 
the energy of the 3-1-2 dissociation channel is higher than that of the 4-1-1 one. 



TABLE III: Energy E of selected exotic five particle systems (in atomic units). The energy E^^ 
of the lowest threshold (low.thr.) is also included, (mp = 1836. 1527 rrj-e, = 3670.4827 me, 
mt = 5496.92158 me, m^ = 206.76826 me) 
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— 1.164023 


lYIIDl 
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— 1 271788 


{ 1 P~ 




—1 149679 


IVIIDl 


(p,e+,e+. 


e",e-) 


-0.8099127 


(p,e+,e- 


,e-) 


-0.788865 


IVllDI 


(d+,e+,e^ 


■,e-,e-) 


-0.81007844 


(d,e+,e- 
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-0.7890280 


IVll El 


(p,p,e-,e 
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-102.750286 


(p,P,At", 
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-102.723336 


lYIIEl 


(d+,t+,^- 


e~, e~) 


-111.889612 


(d+t+,/i- 


-,e-) 


-111.864106 


IVIK^I 


(p,P,P,e" 


,e-) 


-481.605173 


(p,P,P,e 


-) 


-481.580324 



C. (p,p,e ,e ,x+) 

The difference between this case and the previously discussed Hj|'-like systems is that one 
of the heavy particle is different from the other two. Figure ISEl shows how the total energy 
varies with a = mf^/irix). The total energy rapidly decreases toward the energy of the H2 
threshold. The system becomes unbound around maj/mg = 2.5. This result shows that an 
H2 molecule can bind a positively charged particle provided that it is at least 2.5 times 
heavier than an electron. So while the H2 cannot bind a positron it forms a bound system 
with a positive muon /i+ (see Table ITTT|l . The properties of this system are very similar to 
those of Hj]": The two protons and the muon form an isosceles triangle where the two protons 
are somewhat closer to each other than to the muon and the electrons are slightly closer to 
the protons than to the muon. If the mass of approaches the mass of the electron the 
distance between the protons and x'^ increases and the electrons become more and more 
localized around the protons. Eventually, for m^jnif. > 52.5, the system dissociates into H2 
plus x~^. 
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FIG. 26: Energy of (p, p, e , e , x^) as a function of nie/mx- The dots shows the mass ratios where 
the energies were calculated. The horizontal line marks the H2 threshold. Atomic units are used. 

D. {M+,x+,x+,e-,e-) 

By adding two positive charges to the H~ ion, one obtains a x"*", x"*", e~, e~)-like 

system. The stability of this system can be predicted in the limit of very large rrij; by 
noting that the H~-like (M+,e~e~) negative ion behaves like a negative point charge and 
can bind two positive charges, forming a H^-like system. In the limit where rux is very 
small, (M+,e~e~) binds two identical light particles in the same manner as a proton binds 
two electrons in H^. 

The (M"*", x"^, e^, e~) system can dissociate into 4 + 1 [(M+, x+, e^, e^) + x"*" and 
(x+,x+,e-,e-) + M+)] and 3 + 2 [(M+,x+,e-) + (x+,e-) and (x+,x+,e-) + (M+,e-)] 
subsystems. Figure 1771 shows the binding energies as a function of me/m^ < 1, in the case 
where the mass of the heavy particle is equal to the mass of the proton, M = 1836. 1527 me. 
For these configurations with rrie < mx, the lowest threshold is is of the type 4 + 1. 
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FIG. 27: Energy of {M~^,x~^,x^,e , e ) as a function of rrie/rnx (solid line). The dashed line shows 
the energy of the {M^ , ,e~ ,e~) threshold. Atomic units are used. 

Examples of such bound systems are (p, e~, e") or (p, e+, e+, e~, e~), see Table 
mil This latter system will be investigated in detail later. Table IIIll shows that the energy 
of (p, e~, e~), just like that of (p, p, /i"*", e~, e~) in the previous example, is close to 

that of B.^ . The proton and the muon are likely to form an isosceles triangle but now the 
like particles are further away from each other so the base of the triangle is longer than the 
sides in this case. The most important difference is that so long as the rrie/mx ratio lies 
between and 1, this system remains bound. 

E. (M+,M+,e-,e-,x-) 

Another Coulombic five-body system which has attracted attention is the ion. This 
ion is not bound, but the H— H~ potential energy curve has an attractive part beyond 3.5 a.u. 
This leads to speculation about the possibility of resonant states of this system. The fact 
that the is not bound is a consequence of the Pauli principle. Adding a negatively 
charged particle x~ which has the same mass as the electron (but is distinguishable from 
it) to the hydrogen molecule gives a bound system. Its binding energy is about 0.096 a.u. 
The x~ particle can attach itself to the H2 molecule because the Pauli principle does not 
constrain its motion. 

Figure EHl shows the dependence of the binding energy on the mass ratio mx/rrie. The 
threshold in this case is the energy of the (M"*", M"*", e~ ,x~) four-body system. The calcula- 
tion is not trivial because the energy and structure very strongly depend on the m^/M mass 
ratio. For m^jM ^ we practically have a hydrogen molecule. In the case of m^jM ^ 1, 
the (M+,M+,x~) system forms a Ps"^ ion-like system. 
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FIG. 28: Binding energy of (p,p,e~,e~,x~) as a function of rrie/mx. Atomic units are used. 

Due to the heavy masses, the size of this system will be very small compared to that of 
Ps^ and this small (M"*", M"*", x~) system will act as a positive charge and bind the electron. 
The distances between the particles in (M^,M^,x~) will be very small compared to the 
distance between the centre of mass of (M"*", M"*", x~) and the the electron. This system 
can bind one more electron forming (M"*", M"*", x", e~, e~), which is akin to H~. A possible 
choice of x~ is The (p, p, fi~) system is bound, and as the present calculation shows the 
(p, p, e"*", and the (p, p, e~, e~, systems are also bound. These systems remain bound 
even if the masses of the heavy particles are slightly different, e.g., the {M^ , M^, e~, e~, x~) 
system is bound as a rough estimate for 1/3 < M1/M2 < 1. 

F. {M+,M-,m+,m~,x+) 

The next system considered is {M~^ , M~ ,m~^ ,m~ ,171^+). The four-body system 
{M~^ , M~ ,m~^ ,m~) is akin to the hydrogen-antihydrogen system and it is not found to 
be bound if the mass ratio m/M is smaller than about 0.45, see Sec. IVIl If the mass ratio 
m/M is small, the two heavy particle with opposite charges form a small neutral particle 
and the ion formed by the m"*", m~and particle will not be able to form a bound five-body 
system with it. 

G. {M+,M+,M-,m-,m-) 

This system can be characterized by a single mass ratio a = m/M. \{ m < M then 
the dissociation threshold is the energy of the (M"*", M"*", M~, m^) system. The energy of 
{M'',M'',M-,m-) as a function of a is shown in Fig. EHl The {M"^ , M+ , M' ,m- ,m-) 
system is bound with respect to this threshold (see Fig. IHIH) . (p, p, p, e~, e~) is an example 
of this system (see Table UTTj). This shows that a hydrogen molecule is capable of binding an 
antiproton to form a system similar to H~. If m > M then the relevant dissociation threshold 
is given by the energy of (M"*", M"*", m~, m~). The (M"^, M"*", M~ ,m~ ,m~) system is bound 
in the interval 1 < m/M < 2 (see Fig. EH). 
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FIG. 29: Energy of (M+,M+,M ,m ) as a function of m/M for m < M. Atomic units are used 
and it is assumed that m = nie, otherwise the energy unit should be multiphed by m/irie. 
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FIG. 30: Binding energy of (M+,M+,M ,m ,m ) as a function of m/M for m < M. Atomic 
units are used and it is assumed that m = m^ assumed. 
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FIG. 31: Energy of (M'^ , , M~ ,m~ ,m') as a function of m/M for m > M (solid line). The 
dashed line shows the energy of the {M~^ , M'^ ,m~ ,m~) threshold. Atomic units are used and it 
is assumed that m = rriQ. 

There is a very interesting difference between these two cases. In the first case a is 
between and 1. For small a values the three heavy M particles form a positive ion. Its 
size is small and it behaves like a single structureless positively-charged particle . This 
makes it able to bind the two lighter charges forming (c"^, m~, m~). The resulting system is 
very similar to H^. 

In the second case l/cx varies between and 1. Here in the limiting case where a is 
infinite, one has two heavy m~ particles and a composite positive charge "C*"*"" formed 
by (M"*", M"*", M^). This composite particle, however, cannot be viewed as structureless 
in the presence of the heavier m~ particles. Energetically it is more favorable to form a 
(M"^, m~) -|- (M"^, m~) molecule than a (C"*", m~, m~) system so the binding is lost at some 
(T > 1. 

H. (e+, e"^, e^, e^, x^) 

The previous examples started from systems with two heavy positive and two light neg- 
ative charges. The other end of the mass spectrum where one has two light positive and 
two light negative charges has also been investigated. In this case the two negative particles 
were electrons and the two positive particles were positrons. The sign of the charge of rrix 
is not important in this case. For the calculations reported in this section the extra charged 
particle is assumed to be distinguishable from the electron and the positron. 

The five-body binding energy versus the rux/me ratio is shown in Fig. 1221 When the 
system has a mass ratio satisfying rux > rrie, the lowest energy threshold is the energy of 
the (x"*", e~, e~, e"*") -t-e"*") dissociation channel. This system is bound for all mass ratios such 
that rrix > rrie and it can be seen that that the binding energy of the five-particle system 
increases with increasing rrix- 
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FIG. 32: Energy of (m^, m^, m^, m^, x^) (solid line). The horizontal line shows the energy of the 
Ps2 atom, the dashed line is the energy of the (x^ , m'^ , m~ , m~) threshold. Atomic units are used 
and it is assumed that m = rrie. 

When the mass of the distinguishable particle is lighter than that of the electron, i.e., 
rrix < fUe, the threshold energy is the energy of Ps2 + x~^. The binding energy decreases 
steadily as rux is decreased. The system is no longer capable of forming a 5-particle bound 
state when rrix = 0.56 me. The structure of the (e"*", e"*", e~, e~, x"*") system increasingly 
resembles the structure of a system best described as + Ps2. The dissociation limit is 
approached as rrix 0.56 rrie. 

I. The e+PsH system 

The e"''PsH system, (p, e~, e~, e+, e+), corresponds to the case where rrix = rrip, and it is 
clear from Fig. EH that this system is bound. The system is stable against dissociation into 
the H -|- Ps^, p + Ps2 or the PsH -|- e"*" channels. The lowest threshold is the energy of the 
PsH -|- e"*" channel (0. 789197a. u.) and e^PsH is bound by 0.021050 a.u. with respect to this 
threshold (see Fig. [HH]) A model for this system is a positron orbiting the PsH subsystem at 
a relatively large distance from the nucleus. 
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FIG. 33: Energy levels of the HPse"'" and the HPs -|- e"*" and H -|- Ps"^ dissociation channels. 



VIII. LARGER SYSTEMS WITH BOSONS OR FERMIONS 



Molecules are built up of multielectron atoms by various binding mechanisms. The charac- 
teristic property of molecules is that multiple charged heavy atomic nuclei share a delocalised 
electron cloud. The question considered in this section is the existence of bound states of 
a system made of m particles of unit positive charge and n of unit negative charge. Such 
states probably exist in the case of bosons with \m — n\ zero or small. We have investigated 
various systems with m electrons and n positrons, but we were unable to find bound states 
for 4 < m + n < 9. Another possibility that we have investigated is a system comprising a 
proton, m electrons and n positrons, similar to HPs and HPse"*". 

The HPse"*" = (H~, e"*", e"*") is a positively charged system so one may try to add one 
more electron to see if it remains stable. The convergence of the energy is shown in Fig. IH^ 
The energy of the system slowly converges to the lowest (HPs + Ps) threshold and the 
size of the system continuously increases showing that this system is unlikely to be bound. 
Surprisingly, however, by adding two electrons to (H~,e"'",e"'") one gets a bound system, as 
shown in Fig. OH This system "H~Ps2" contains a proton, two positrons and four electrons, 
and can also be considered as a three-body system made up of a proton and two Ps~ ions 
analogous to the H~ ion (with the electrons replaced by composite Ps~ ions). 
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FIG. 34: Convergence of the energy of the (p+,e+,e+,e^,e^,e^) and (p+,e+,e"'",e~,e~,e~,e^) sys- 
tems. The dotted line is the HPs + Ps, the dashed Une is the H~ -|- Ps2, the long dashed line is the 
HPs + Ps" threshold. 

In the case of neutral, self conjugate, bosonic system s (rn ^)", (m")", where the mass can 
be chosen as m = 1, a crude variational calculation |l83j gives energies -E(3) = —0 .789 
and E{A) = -1.046. An improved value E{3) = -0.820 was obtained in Ref. |l8^ A 
large-n behaviour E{n) oc n^/^ has been established |l85j . and heu ristic arguments jl86j 
suggest that E{n) ^ — 0.148 n''/^. The exact or approximate |l83l | values correspond to 
E{n)/nJ/^ = —0.25, —0.196, —0.169, —0.150 for n = 1, ... ,4, respectively. It is not known 
whether E{n)/n'/^ is expected to be converge monotonically to the value suggested by 
heuristic arguments. If so, it is remarkable that it is so close to the suggested limiting value 
for n as low as 4. 



IX. SYSTEMS CONTAINING ANTIPROTONS 



This review mostly deals with how small atoms and molecules evolve when the familiar 
negatively-charged electrons and positively-charged nuclei are replaced by other elementary 
constituents. Investigations with positrons, and to a less extent, pions, muons or kaons 
have been carried out for many years, and are in the process of further development. In 
recent years, however, the most spectacular progress have been observed in the domain of 
antiprotons, thanks to the commissioning of sources of intense and cooled antiprotons, in 
particular at the European Laboratory for Particle Physics (CERN). In this section, we shall 
review the current state of knowledge in this area and prospects for the future. 



A. Antihydrogen, protonium and other antiprotonic atoms 

1. Antihydrogen atoms 
Antihydrogen, H = (e"*", p), is the charge conjugat e of the familiar hydrogen atoms. It has 



been produced for the first time at CERN |l87l . Il8q . by directing a high-energy antiproton 



beam onto the electric field of a nucleus, where it enc ountered virtual (e"'",e~) pairs. The 
experiment was successfully repeated at Fermilab [l89'|. 

Antihydrogen production is of considerable current interest on account of the successful 
preparation at CERN in 2002 of about 50 000 atoms of H atoms, cooled to within 15 
degrees of absolute zero. See Amoretti et al. [190]. This followed on from several years 
preparatory work by the ATHENA project (ApparaTus for High pre cisio n Experiments 
on Neutral Antimatter), Charlton et al. 1 191] and Holzscheiter et al. |l92l Il93j . and the 
competing ATRAP (Antihydrogen TRAP) project of Gabrielse et al. The most recent 

results obtained by ATRAP are described by Gabrielse et al. |l95j . 



Now that cold H atoms have been prepared, it should be possible to trap H at a sufficiently 
low temperature that the laboratory frame is essentially is essentially its rest frame. This 
would make possible tests of the predictions of two fundamental theories of modern physics: 
quantum field theory and Einstein's general theory of relativity. In particular, tests could 
be made of the charge, parity and ti me-re versal (CPT or TCP) symmetry of quantum 
field theory. See, for example, Hughes |l96| . Schweber jl97l l. and for a recent discussion, 
Ref. (l98| . 

It follows from the CPT symmetry that a charged particle and its antiparticle should 
have equal and opposite charges. Also their masses and gyromagnetic ratios should be equal 
and their lifetimes should be the same. In additi on, this symmetry predicts that H and H 
should have identical spectra |l9lL ll92L Il93l . Il94| . Experimentalists plan to test, as far as 
possible, whether H and H do have these properties. In particular, they intend to compare 
the frequency of the ls-2s two-photon transition in H and H. 

The equivalence principle, according to which all bodies fall at the same rate in a gravi- 
tational field, led Einstein to his general theory of relativity in which gravitation manifests 
itself as a metric effect, curved spacetime, as opposed to the flat spacetime of special rel- 
ativity in the absence of gravitation. See, for example, Weinberg [199] and Hughes jl96j. 
The exact form of the curvature is determined by Einstein's fleld equations. 

It is proposed to test the validity of the equivalence principle by carrying out a null red 
shift experiment in which the frequency of the two-photon ls-2s transition is observed for H 
and H as both are moved through the same gravitational field. Any difference between the 



two sets of values would indicate that the two atoms experienced a different gravitational red 



shift, which can be shown to be a violation of the equivalence principle 191 . 196l |. However, it 



must be said that the possibility of such a violation is much more remote when gravitation is 
viewed from the standpoint of general relativity as resulting from the curvature of spacetime, 
rather than as a force obeying an attractive inverse square law, as in the Newtonian theory. 



2. Protonium 

Protonium, (p, p) often denoted by Pn, is the heavy-mass analogue of the positronium 
atom. The Is energy and radius are of the order of 12 keV and 50 m. 

When antiprotons are stopped in a dilute hydrogen target, antiproton capture takes 
place in a state of high principal quantum number n of the order of n = 40. The electron 
is expelled and the antiproton quickly cascades down to lower orbits, with circular orbit 
i = n — 1 preferentially populated. Strong interaction effects are negligible as long as n > 2. 
For the 2p state, strong interaction becomes important. States are shifted, but due to the 
pronounced spin dependence of the nucleon-antinucleon interaction, the singlet state, and 
the triplet states with J = 2, 1 or are moved differently. The annihilation with also 
becomes dominant with respect to radiation to the Is ground-state. If the target is chosen 
to be more dense, protonium states often penetrate inside the hydrogen atoms, where they 
experience a strong electric field which produces a Stark mixing of {n, i = n — 1) states with 
states of lower i, and eventually states with £ = or £ = 1, in which annihilation takes 
place. This is the Day-Sucher-Snow effect, which is important to understand the initial 
state of a nnih ilation experiment. For a recent review on protonium, and references, see, 
e.g., Ref. [iooj . 

Exotic atoms such as protonium give rise to an interesting quantum-mechanical problem, 
in which the Coulomb potential —Z/r is supplemented by a short-range term U{r), which is 
sometimes chosen to be complex, to mimic the effect of absorption or annihilation. A naive 
application of perturbation theory suggests that the energy E changes with respect to the 
pure Coulomb energy Ec (with wave function \E'c) by 

E-Ec ~ y |^(°)p?7(r)dV , (9.1) 

which is unsatisfactory, as it gives too large a correction. The shift, indeed, remains small, 
even if U{r) is an infinite hard core. The shift, in fact, is proportional not to the strong- 
interaction potential U, but to its scattering length a, and this shift is small as long as |a| 
remains small compared to the Bohr radius R. More precisely for the n-s statc[200], 

E-E,^ -E,- ^ . (9.2) 



3. Ordinary antiprotonic atoms 

When an antiproton is captured in the field of a nucleus, it usually cascades down very 
quickly toward low-lying orbits close to the nucleus. The electrons surrounding the nucleus 
are either expelled or remain at large distances. Hence, with the remarkable exception 
of metastable antiprotonic Helium states, on which more below, the few-charge problem 



decouples: within a distance of a few fermis, or a few tens of fermis, we have a two-body 
atom consisting of a nucleus of charge Z and an antiproton; then, at the usual atomic scale, 
we have an ion with an effective nucleus of charge {Z — 1) and a few electrons. 

The energy levels of the antiproton are shifted and broadened by strong interaction. 
This gives a measure of the antiproton-nucleus interaction, to be compared with data from 
antiproton-nucleus scattering experiment, and theoretical expectations based on folding the 
antiproton-nucleon amplitude w ith t he wave-function of the nucleus. For a review and 



references, see, for instance, Ref. |201 



B. Three unit-charge systems with antiprotons 

If one combines and antiproton with protons, positrons, and electrons, one gets the 
following three-body systems. 

• (p, e"*", e"*") is the conjugate of H^. It could in principle be formed and compared to 

to test CPT symmetry, 

• (p,e"'",e~) is unstable, 

• (P) P; p) is stable with Coulomb forces only. Its study could probe the long-range part 
of the strong interaction between a proton and and antiproton. 



C. Antiprotonic helium 

In the limit where the positive charges are very heavy compared to the negative 
charges, the four-charge configurations {M^ , ,m~ ,m'~) reduce to Helium-like atoms, 
say {M~^^ ,m~ ,m'~). The spectroscopy of the He atom is well known and documented 
|l6|, and thus will not be discussed further. The existence of at least one bound state of 
(M^"*", m~, m'~) is rather intuitive: the M~^^ nucleus can bind the heaviest of m and m', 
and there is a long-range Coulomb force left to bind the lightest negative charge. 

Let us take first the limit where M is infinite. If m = m', we have the familiar Helium 
atom, with the two electrons in the Is orbit. Indeed, unlike the case of H, for which the 
correlation between the electrons is crucial (see Sees. IIII Hl and lD^ . If, for instance, m' 3> m, 
the a particle and the particle of mass m' form a compact core of charge which forms 
a hydrogen-like atom with the light particle m. The binding of m and that of m' occur at 
different scales, and there is almost no three-body effect for this ground state. 

Kondo, and Russel, inspired by anomalies in early data on kaonic atoms, have shown that 
the dynamics becomes more intricate for some excited states. For references, see, e.g., j20^. 
If the heavy negative charge m' is in an excited state, such that its mean separation from 
M^"*" becomes comparable to that of m, then the radiative and Auger transitions become 
suppressed. Experiments with negative muons and later with antiprotons on Helium 
have revealed, indeed, delayed transitions to lower states. 

The utility of these long-lived states is even better than might be expected. Some of 
these antiprotonic "atomcules" have smaller intrinsic width than conventional excited states 
of ordinary atoms, and thus can serve as a improved prototypes to define and measure 
fundamental quantities (antiproton mass or magnetic moment, fine structure constant, etc.). 



It is not our aim here to review all properties of these antiprotonic Helium states, nor 
to discuss the elaborate theoretical calculations by Ko robo v et al. and Y.Kino et al. We 
refer the reader to t he rec ent review by Yamazaki et al. 202| |. and to the Proceedings of the 
LEAP03 conference ET" 



In the spirit of our discussions on the how the stability of three- or four-charge systems 
depends on the involved masses, we note that the metastability of these new atomcules 
has been studied with ^He and ^He, and also with muons and antiprotons.The use of an- 
tideuterons j204| . The CERN antiproton source could be tuned to produce, collect and 
cool antideuterons, at a rate of about 10^'^ that of antiprotons, largely sufficient for capture 
experiments. 

An interesting question is also whether the dramatic metastability observed in CERN 
experiments could survive a smearing of the Helium core. Suppose one replace He"*"^ by two 
deuterons: are there metastable states in the four-body molecule (d, d, p,e~)? 



D. Four unit-charge systems with antiprotons 



1. Symmetry tests 



The antihydrogen molecule (p, p, e"*", e+), the protonium molecule (p,p,p,p) and the 
positronium antihydride (p, e~, e"*", e"*") are examples of configurations that are stable if the 
interaction is purely Coulombic. They will only become interesting if they ca n be used in 
high-precision experiments to test the CPT symmetry of quantum field theory |205j | or mea- 
sure the proton-antiproton interaction at a range of distances at which it cannot be obtained 
accurately from a spectroscopic study of protonium. 



2. HH 



There has been in teres t in calculations on t he in teraction of H with H, since the 1970s, 
Junker and Bardsley |206|. M organ and Hughes |207j , Kolos et al. |208j , Campeanu and Beu 
|209j , Shylapnikov et al. |21[| . However, the work on the preparation of H described above 
has gene rated greater theoretica l interest i n th i s system. Armour al. |l69L 211 . bid l21," ' 



irim^j Jonsell, Froelich et al. ^^^^^^^^^M,^^^^, Sinha and Ghosh ^ 
Voronin and Carbonell |224| . 

The H-H system is a four-body system similar to the H2 molecule, in that it contains 
two nuclei of the same mass as the proton and two light particles with the same mass as the 
electron. However, both the nuclei and light particles have equal but opposite charges. As 
might be expected, this leads to very significant differences from H2. 

One of the most obvious is that if the nuclei coincide in H2, the resulting system is He, 
whereas when this happens in H-H, the resulting nucleus has zero overall charge and thus 
clearly cannot bind the light particles. 

It is reasonable to assume that the critical internuclear distance Rc, below which a fixed 
proton and antiproton cannot bind the electron and the positron is greater than 0.639ao. 
This can be seen as follows. We know from the work on (p, p, e~) described in Sec. IIIIFI and 
the invariance of the energy levels of this system if the electron is replaced by a positron 
that, if i? < i?c; the electron and the positron would be unbound if they did not attract 



each other. In addition, the attraction between them makes possible separation into ground 
state positronium, which is — 0.25a.u. lower in energy. 

It might be thought that such a result would be easy to prove. However, this does not 
seem to be the case. 

Kolos et al. [208] showed that Rc < 0.95 ao, using the Ravleigh-Rit z va riational method 
and a trial function containing Hylleraas-type functions. Armour et al. |212l | went on to show 
R < 0.8 ao, a fixed proton and an antiproton can bind an electron and a positron. They used 
a method similar to Kolos et al., but included a basis function to represent weakly bound Ps, 
a long way from the nuclei. More recently, Strasburger |.225i] has shown that Rc < 0.744 ao, 
using the same variational method but with a basis set made up of explicitly correlated 
Gaussian functions, similar in type to those used by Kamimura 117| in his calculation of 
the energy levels of (d,/i~,t). 

The lowest continuum threshold of H-H is Pn + Ps, where both are at rest in their ground 
states, at infinite separation. Note that we can obtain continuum states with any angular 
momentum, as close in energy as we please to this threshold, by associating the angular 
momentum with the motion of the centre of mass of either the Pn or the Ps. The threshold 
has a very low energy, — 459. 29 a. u. This is due to the very low energy of the ground state 
of Pn. For Pn in this state, the expectation value of R is 0.0016 ao, which is far below the 
value of Rc. In view of these two properties of Pn + Ps, it is very u nlike ly that a bound 
state of H-H exists, though this has not so far been proved rigorously |l67l |. 

Mo re d efinite evidence for the absence of an H-H state has been obtained by Bressanini 
et al. jl52| using the variational Monte Carlo method (VMC) |22fij . They investigated the 
stability of systems made up of two pairs of particles (M+, M~) and (m"*", m"*"), where each 
pair has the same mass and charge ±e, where e is the charge on the proton. According 
to the VMC method, if M is chosen so that M > m, then no bound state of the system 
exists if M/m > 1.7. Using the diffusion quantum Monte Carlo (DMC) method, they obtain 
M/m > 2.2. These results are likely to be accurate but are not rigorous. However, in view 
of the fact that M/m = 1836 for H-H, it seems incredible that there should be a bound 
state in this case. The high value of the mass ratio strongly suggests that H-H must lie far 
inside the domain of instability of systems of the above type. 



E. Antihydrogen Helium 

We have seen that in the case of HH if the nuclei are fixed then there is a critical value 
Rc, of the internuclear distance, R, of the order of a Bohr radius, below which the nuclei 
are u nable to bind the electron and the positron. The situation is different in the case of 
HeH |214| . This is because the He nucleus has charge +2. Thus the nuclear potential is 
Coulombic asymptotically, rather than dipolar as in the case of HH. It is attractive for the 
electron and repulsive for the positron. 

At i? = 0, HeH corresponds to positro nium hydride (PsH). This is known to have a 
bound state with binding energy, 1.06 eV j227l |. There are continuum thresholds for HeH 
corresponding to antiprotonic helium + positronium, He'''p+Ps, and Hep -|- a positron (e"*"), 
as well other thresholds involving more extensive decomposition of HeH. 

At i? = 0, He'''p corresponds to H + Ps and Hep + a positron to H~+ a positron. They 
have energies of —20.41 eV and —14.45 eV, respectively, with respect to completely separated 
systems. Thus the the continuum threshold involving antiprotonic helium is lower at i? = 0. 



However, as R 



He+p + Ps ^ He+ + p + Ps , (9.3) 

whereas 

Hep + e+ ^ He + p + e+ . (9.4) 

In this hmit, the continuum threshold associated with the antiprotonic hehum plus positro- 
nium has an energy of — 61.2eV, whereas the threshold associated with Hep has an energy 
of -79.0 eV. 

It follows that these two thresholds must cross at some point between R = and R = oo. 
As the Schrodinger equation for the Born-Oppenheimer potential of He'''p helium is separable 
in prolate s pher oidal coordinates, the threshold energy for He'''p + Ps can be calculated very 
accurately l228|. The Born-Oppenheimer potential for Hep can be calculated accurately us- 
ing the Rayleigh-Ritz variational method. The results obtained for both thresholds indicate 
that crossing occurs in the vicinity of R = 1.25 ao |22^. A ve ry ac curate Born-Oppenheimer 
potential has been obtained by Strasburger and Chojnacki 23Cl| |. using the Rayleigh-Ritz 



variational method and a basis set made up of 768 explicitly correlated Gaussian functions. 
Comparison of the value of this potential with the energies of the two continuum thresholds 
considered above, shows that the helium nucleus (a) and the antiproton are able to bind the 
electrons and the positron at all R values. This is in contrast with the situation in the case 
of HH in which the electron and the positron become unbound if R is less than the critical 
value. Re- 
in fact, of course, the nuclei are not fixed. We have seen that in the case HH the ground 
state energy of protonium is very low and the expectation value of R in this state is very 
small, so small that is incredible that the dipole formed by the proton and the antiproton 
would be able to bind the electron and the positron. 

The ground state of (a,p) has an energy of — 2938. 2 a. u. and the expectation value of 
R in this state is 0.0005 ao, both of which are considerably lower than the corresponding 
values for HH. However, as noted earlier, there is a key difference between the two systems: 
whereas Pn forms a dipole, a + p has a net charge of +1. Thus, it is very close to being a 
single particle with the same charge as the proton but a much larger mass. 

As mentioned earlier, H + Ps is known to form a bound state, positronium hydride (PsH). 
If the inclusion of the three leptons only perturbs the ground state of the nuclei to a small 
extent, so that there is a state of Hep in which the a-nucleus -|- p are in a state that differs 
only slightly from the ground state of these nuclei on their own, then it can be expected 
that HeH will have a ground state whose leptonic wave function is very similar to that of 
PsH. It would thus have a bound state. Furthermore, removal of a positron from the system 
in this state would give rise to Hep in a state with electronic wave function very similar to 
H~ in its ground state, which we have seen earlier is a bound state. Thus, if the presence 
of the leptons has only a small effect on the nuclear ground state, both HeH and Hep will 
form bound states. It will be of interest to see whether the existence of such states can be 
proved. For a recent study of HeH, see |23ll | . 



F. Perspectives 

Positrons are nowadays commonly used in fundamental and applied physics. The physics 
of antiprotons is becoming accessible thanks to the development of sophisticated trapping 
and cooling devices. 



Light antinuclei (d, t, ^H e) ar e currently produced in antiprotons factories (d was identified 
in 1965 by Zichichi et al. (232| |). The rate is much lower than for p, but is is sufficient to 
form and study a few atoms. 

So, in principle, one could conceive of delicate systems with d instead of p or a mixture 
of p and d, such as the Borromean molecule discussed in Sec. I VI (^1 



X. FEW-CHARGE SYSTEMS IN TWO DIMENSIONS 



Two-dimensional (2D) charged-particle systems are often used as a model of excitonic 
complexes confined in semiconductor quantum wells. When a thin layer of a material on the 
xy plane is sandwiched between two layers of another material that has a different dielectric 
constant, electrons and holes in the middle layer are basically confined by a potential Vc{z) = 
if \z\ < d/2 and Vc{z) = Vq otherwise, where d is the thickness of the middle layer, with a 
large value of the confining potential Vq. Hence every wave-function approximately factorises 
into a z-component which is the ground state of Vc, and an horizontal component where the 
interesting dynamics takes place. This is why 2D-systems are thoroughly studied, in addition 
to their intrinsic interest. Some basic properties of bound states in 2D-quantum mechanics 
are reminded in Appendix I A 101 

The most important property of the confined 2D system is the increased binding en- 
ergy. This is even more pronounced for large systems than for the small ones of which 
the dissociation threshold consists, and the binding energy increases as compared to the 
three-dimensional case, even if measured in units of the threshold energy. 

The binding energies of some 2D and 3D electron-hole systems are compared in Table IIVI 

TABLE IV: Energies and binding energies of the 2D and 3D excitonic complexes in the hydrogenic 
(fj = 0) and positronium (a = 1) limits, a = m*/m'^ is the ratio of electron and hole effective 
masses. An asterisk indicates that the state is fomid to be unbound. The unit of energy is the 
excitonic Rydberg. 
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The total energies of 2D electron- hole systems as a function of the mass ratio a = me/mh_ 
are shown in Fig. 123 (see Fig. |2Slfor the corresponding 3D case). The binding energy of the 
2D system is almost by a factor of 10 larger than that of the corresponding 3D system. 




FIG. 35: The total energies of 2D electron- hole systems as a function of the electron hole mass 
ratio a = me/m\i. 

The binding energy in general decreases by increasing a from to 1, and this trend is 
especially dramatic in the case of X^. The 2D {eeh), [ehh) and (eehh) (Fig. refvkf9 systems 
are bound irrespective of the mass ratio, just as in the corresponding 3D cases. In 2D and 
3D the (eehhh) forms a bound system for small values of cr, while the (eeehh) system is 
unbound for any values of cr. The critical value of a is 

cTcr < 0.27 for 2D , , , 

10.1 

cTcr < 0.23 for 2D . ^ ^ 

The properties of the 2D and 3D systems are found to be generally very similar. Even 
though the binding energy and therefore the relative distances between particles are very 
different in 2D and 3D, the binding energies show very similar behaviour as a function of a 
(see e.g.. Fig. It is a striking similarity that the stability of the (e, e, h, h, h) is lost at 
nearly the same value of acr in 2D and 3D. 

We should point out here that the electron-hole complexes are highly nontrivial quantum 
mechanical systems in both 2D and 3D. The relative motion of the particles is complicated 
and it is impossible to model these systems using some rigid geometrical picture, e.g., by 
assuming that the (e, e, h, h) forms a static square in 2D. The interpretation of the (e, e, h, h) 
molecule as a system of two exciton atoms is also an oversimplification. 




FIG. 36: The binding energy of the biexciton X2 (eehh) compared to the binding energy of the 
exciton X as a fu nctio n of the mass ratio a. The dashed curves refer to a variational calculation 
described in Ref. |233i |. 



XI. CONCLUSION AND OUTLOOK 



In this review, we have addressed the question whether a system of point-hke charges, 
interacting through their Coulomb potential, form a collective bound state with all particles 
remaining at finite distance the ones from the others, or split into clusters or isolated charges. 

In a few cases, the answer is obvious: a positive charge q = +2 can bind two electrons, 
since after having attached a first electron, it still exerts a Coulomb attraction on the second 
one. On the other hand, a tiny positive charge g — > cannot bind two electrons. 

In most case, the answer depends on the mass ratios. For three unit charges, the positro- 
nium ion (e"'",e~,e~) and other symmetric systems are stable, as well as nearly symmetric 
systems, while very asymmetric systems such as (e~, p, e"*") or (p, p, e^) do not form a three- 
body bound state. 

For the case of four unit charges, {m^ , m^, ,171^), there is also a premium to symmetry, 
since all systems with two identical particles, e.g., ms = m4, are stable, and among them 
the twice-symmetric positronium molecule, Ps2 = (e"*", e"*", e~, e~). However, starting from 
Ps2, it is observed that breaking charge-conjugation symmetry into (M"*", M"*", m~ ,m~) with 
M improves stability, whilst breaking permutation symmetry into (M"*", m"*", M~, m~) 
quickly result into a loss of stability. The question is whether symmetry breaking benefits 
more to the collective state than to its threshold, or vice-versa. 

This problem of stability has been investigated by many authors since the beginning of 
quantum theory. Rigorous results have been obtained, with simple or sophisticated methods 
of mathematical physics. This subject has also stimulated extremely accurate and innovative 
variational calculations. 

In nuclear physics, there are striking differences between compact, double-magic nuclei, 
where nucleons are tightly bound together, and halo nuclei with a tail of neutrons extending 
very far, though in both cases, we are dealing with stable bound states. One can similarly 
notice wide differences among bound states of charged particles. In the case of three-body 
bound states with unit charges, it is customary to distinguish "molecular" from "atomic" 
ions. The former ones, whose prototype is {m^ , m2 , m'^) = (e~,p,p) = H^, ha ve a large 
excess of binding below the threshold, easily survive symmetry breaking m2 7^ rn^ |234| , and 
have several excited states. The latter ones, as the hydrogen ion H~ = (p, e~,e~) have not 
stable excited state, and a small binding energy which quickly disappear if tt t,? is allowed 
to differ from 7713. The change from aton aic t o molecular states is analysed in |l2fil |. with a 
different language, and also identified in |234| . 

For four-body systems, there is even a wider variety of situations. While Ps2 is weakly 
bound, the H2 ground-state lies more deeply below its threshold. Some states are well 
described as two atoms with a relative orbital motion, while others look more as a compact 
ion surrounded by a light particle. After the completion of this review, an investigation of 
(m"*"^, e+, e~, e~) stability was investigated by Mitroy and Novikov |235j | . with the mass and 
charge of first particle being varied. 

Clearly several aspects deserve further investigations, in particular for systems of more 
than four particles, for which the available numerical investigations could be extended and 
summarised by outlining the shape of the stability domain. 

Many of the results presented here can be extended to cases where the potential is not of 
Coulomb type. Let us give an example. When discussing the shape of the stability domain 
for a 3-body system {mf ,m^), we have used the property that in the sub-domain with 
the lowest threshold (m^, m^) stability survives and is even improved when the mass 



m3 is increased while mi and m2 are kept fixed. This result holds for any 3-body Hamiltonian 
with a potential that does not depend on the constituent masses. 

The Coulomb and other power-law (r'^) potentials provide the simplification that the 
energies scale like E oc m~^/'^^^'^\ when a factor m is applied to the masses. Hence the 
problem of stability of a A^-body system depends only on [N — 1) mass ratios. 

There is, however, a major difference between on one hand Coulomb or other long-range 
potentials and on the other hand dipole-dipole potentials or sharply decreasing potentials 
such as the Yukawa potential: the former support a countable infinity of bound state however 
small is the coupling, while the latter require a minimal strength to achieve binding. Note 
that, as pointed out in Sec. IIII Fl once the minimal strength has been attained, a dipole- 
dipole potential can support a countable infinity of of bound states. 

The transition from short-range to long-range potentials is interesting. Bressanini et al. 

|286j . for instance, studied how the spectrum of familiar atoms and molecules evolveo when 
the Coulomb potential becomes screened, say 1/r ^ exp(— Ar)/r. Interestingly, they found 
that larger systems are better armed to survive screening. For some values of the range pa- 
rameter A, the H2 molecule remains bound while all smaller systems have lost their stability. 
This means that the Borromean behaviour denoted for the exotic configuration (p, d, p, d) 
is already underneath for more familiar ions and molecules. The subject of Borromean con- 
figurations with more than four constituents, with or without screening, is one of the many 
developments to be envisaged in this domain. 



APPENDIX A: SOME BASIC RESULTS ON SCHRODINGER OPERATORS 



We summarise here some results that are often referred to in this review. 



1. Variational principle for the ground-state 

The best known formulation reads as follows. For a system with Hamiltonian, H, the 
expectation value of the energy of any trial wave function, $, is an upper bound to the 
ground-state energy, Eq, i.e., if \E'o is the normalised ground-state wave-function, then 

Eo = i^oim^o) < Em ^ ^-j^ . (Al) 

This result if often used to prove binding, as to prove this, it is sufficient to obtain one wave 
function $ such that E[^] lies below the appropriate threshold. It is also the basis of the 
Rayleigh-Ritz variational m etho d, where $ is parametrised to approach \E'o as closely as 
possible. See, e.g., Refs. |i3,E3- 

The variational principle implies the following stationary property, which is valid for all 
levels. If $ = + where cp is small, $ will be "close" to and E[^] will differ from 
En by a term that is second order in the small error function, cp. 

Note that if one separates out the cm. motion, say H = {J^PiY / + -^int, the 

variational principle with square- integrable wave function applies to Hi^t, but, in practice 
1$) = ($|ifint|$), if one deals with a trial function $ which depends only upon relative 
distances. In short, there is no need to remove explicitly the centre-of-mass motion. 



2. Variational principle for excited states 

The variational principle keeps the same simple formulation for the lowest level differ- 
ing from the ground state by any quantum number corresponding to an exact symmetry. 
Examples are the first state with angular momentum J, or the first state with odd parity, 
etc. 

As noted, e.g., by Karl and Novikov |237l | . one can still use the stationary property to 
get reasonable estimates for radial excitations. For instance, a two-parameter trial-function 
with a node would give a good idea of the radially excited s-wave in a central potential. 
However, more parameters would lead to a wave function with an almost invisible node near 
r = that would approach the ground-state. In other words, there is no guarantee that the 
wave function would remain close to the 2s state during the minimisation. 

To get an upper bound and achieve converged variational calculation, one used the "min- 
imax" principle. The n-th eigenvalue (in a sector of given conserved quantum numbers) is 
bound by the largest eigenvalue of the restriction of H to an n-dimensional subspace Ttn of 
the Hilbert space H, 

En < max ^ ', ' ' . (A2) 

Diagonalising the restriction of H to Tin leads to an upper limit of the n first levels. The 
parameters defining Hn can be tuned to approximate best a given level, and then used 



to describe simultaneously several levels. This ensures the orthogonality between the ap- 
proximate wave functions, a requirement which is important for estimating some transition 
rates. 

In practice, Tin is often constructed by adding a new term in the set of ba sis fu nction 



defining TCn-i, and hence TCn-i C Tin- The Hylleraas-Undheim theorem [121, Il44| states 
that the eigenvalues of Ef''^ obtained by diagonalising the restriction of H to Tin intertwine 
the E^^~^^ obtained from Hn-i, 

...< ) < < eII\ < . . . (A3) 

A consequence of Eq. flA2|) is that if there are n linearly-independent approximate eigen- 
function with energy expectation below the threshold energy E'th, then the number of bound 
states, iVb, fulfils 

iVb > n . (A4) 



3. Comparison Theorem 

If an upper bound on the number of bound states is required, a lower rather than an 
upper bound to the eigenvalues of the system is required. This is, in general, much more 
difficult to obtain [78^ . 

The comparison theorem (see, for example, Ref. jl^]) is a very useful method for obtaining 
an upper bound on the number of bound states of a system. This theorem states that if two 
Hamiltonians Hi and H2 are such that 

(^l^il^) < (^|i^2|*) , (A5) 
for any allowed square- integrable function \E', then 

E[i! < , (A6) 

where E^^ is the energy at which the continuum threshold (if any) of Hi begins. Also if Hi 
and H2 have at least n bound states and the mth bound state of Hi has energy Em , 

E^^^<Ei^ m<n. (A7) 
In addition if Hi and H2 have the same continuum threshold, i.e., 

E^I! = , (A8) 
and Hi has exactly n bound states, n is an upper bound to the number of bound states of 

In practical applications of this theorem, H2 is taken to be the original Hamiltonian for 
the problem, while Hi is the Hamiltonian of a more tractable system. It follows from this 
theorem that if e[^ = E^.^ , and Hi has no bound states then H2 also has no bound states. 



4. Scaling 



Consider a typical Coulomb Hamiltonian 

^ 2m,- ^ ' ^ ^ 

Kj -> 



and multiply all masses by m and all charges by q, and expand all distances by a factor p, 
in particular rjj —>■ prij. The kinetic term gets a factor m~^p~'^ and the potential ones a 
factor q^p~^. These factors become equal to an overall energy rescaling by a factor e = mg^, 
if the distances are scaled by p = m~^q~'^. 

This scaling can also be understood from the stationary property of the eigenstates. The 
expectation value of H has a vanishing derivative at p = 1, if the exact wave function 
\l/(ri2, . . .) is used to build a tr ial wave function = "^{pr 12, . . .)• 



As noted by Hylleraas [130[, Fock |131l | and many others, most basic properties of the 
exact solutions remain true for variational solutions, with the mild restriction that the set 
of trial wave functions is globally invariant under rescaling. This is true for the scaling 
properties of energies and wave functions. 

Scaling considerably reduces the parameter dependence of the few-body Coulomb prob- 
lem. In particular, all hydrogenic systems 

H-.nt = ^--, (AlO) 

have similar energy spectra and eigenfunctions. The properties of a Coulomb system con- 
taining N particles depend only on (A^ — 1) mass ratios and (A^ — 1) charge ratios. 

Scaling laws are easily generalised to other power-law potentials rf^. See, e.g., Ref. |238l |. 

Partial scaling holds for most other potentials. Consider for instance the spectrum of 

^int = ^ - g'-^^tl^ . (All) 

This Hamiltonian is proportional to 

= - G^2Htl) . G=^. (A12) 

and thus, it is sufficient to know its properties for all but one of the parameters h, p, m,, g 
set to 1. 



5. Virial theorem 

For a bound state in a Coulomb Hamiltonian ()A9jl . the energy E = (H), is the sum of a 
kinetic term T = (^pf / {2mi)) and a potential energy U = (^qiqi/rij) , such that 

T = -| . (A13) 

With the same conditions as for scaling, this property extends to variational estimates, thus 
reducing by one the number of parameters to be varied (see, e.g.. Appendix iDl below) . 

For other power-law potentials, the rule is T = aU /2 [238|. For more general potentials, 
the potential energy all is replaced by the expectation value of rd\^/dr. 



6. Hellmann Feynman theorem 



Suppose H depends on a parameter A and let E{X) be a particular eigenvalue, with 
normalised wave function ^E'(A). Then 

^ = (*(A)|^|^(A)). (AM) 

7. Concavity property of the ground state 

We further assume that H depends on A linearly. Then, the ground-state energy is 
a concave function of A. Indeed, the second order derivative, as given by second-order 
perturbation theory, is always negative. Another, simpler, derivation makes use of the 
variational principle [89^ . 

The result does not hold for individual radial excitations, but it is true for the sum of n 
first levels. 

Again, the result can also be written for variational solutions, provided the same set of 
trial functions is used over the whole range of values of the parameter A. 



8. Symmetry breaking 

We often use the property that the ground-state of an Hamiltonian H with a certain 
symmetry has its energy "lowered" if a symmetry-breaking term is added to H. 

Consider first the one-dimensional oscillator (with h = 1) ho = + x^, whose ground 
state lies at eo = 1. It becomes e(A) = 1 — A^/4 if ho ^ + Ax, as can be seen by direct 
computation. 

More generally, if 

H{\) = i/evcn + Ai/odd , (A15) 

is an even Hamiltonian modified by an odd term, one can use the variational principle with 
the even ground-state ^1/(0) of -ff(O) as trial wave function, and derive the following relation 
for the ground-state energy 



E{\) < (M/(0)|/7(A)|M/(0)) = (vl/(0)|iJeven|^(0)) = i?(0) . 



(A16) 



Examples are given in this review where parity is replaced by charge conjugation or per- 
mutation symmetry. Consider, for instance, an exotic Helium-like atom (a,/i~,7r~), in the 
limit where m^, = oo. A symmetric (a,yU~,yU~) or (a,7r~,7r~) system is easily deduced from 
the familiar benchmark (a,e~,e~), by a scaling factor m^/rrie or m.„/me. The (a,/i~,7r~) 
Hamiltonian can be written as 



if (a, /i 
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2 2 



2 2 1 

+ — + 5 (m~M 

ri r2 ri2 



Pi 



Pi 



(A17) 



where m' 



{m^ + m^^) /2 and 5{m ^) = { 



m 



-1 



m^^)/2, and hence the ground state of 



{a,fi \ TT ^) is close, but slightly below a rescaled (a, e ,e ), i.e. 



TT 



< E{a,e ,e" 



(A18) 



the difference being of second order in 6{m~^). A similar averaging of inverse masses is 
adequate for isospin-breaking effects in simple nuclear models or quark models with isospin- 
conserving forces, to account for the mass difference between proton and neutron or between 
u-quark and d-quark. 

A current but not always correct statement is that "a symmetry-violating term will push 
the lowest (even) level down and the first excited (odd) state up". This is true if the 
additional interaction has negligible matrix elements between any of two first levels and the 
higher states. Then the effect of symmetry breaking can be described by a 2 x 2 matrix, 




(A19) 



and the sum of eigenvalues (trace) is conserved as the mixing strength e is varied. This 
is, however, a closed system, a condition that often does not hold for physical examples. 
For instance, in the case of + + Ax, all levels are pushed down when the odd term is 
switched on. 



(A20) 



9. Jacobi coordinates 

For two-particles systems, the transformation 

{t-i, raj ^ {r = r2 - ri, R= ■ } , 

mi + m2 

P2} ^ {P = : , -P = Pi + P2} , 

nil + 

makes it possible to separate out the motion of the centre of mass, 

— H = \ , mi2 = nil +1712 , ^h_2 = f^i^ + "^2 ^ • (^21) 

nil fn2 nii2 jjL 

For three particles, one can apply this transformation to particles 1 and 2, and then to 
particle 3 and the centre of mass of 1 and 2. This means introducing 

/o = T2 - ri , A oc rg ■ , R = — , (A22) 

nil + ^2 



as illustrated in Fig. EH and the conjugate momenta, satisfying the relation 

vii. ii + f!i + ^. (A23) 



The coefficient of A can be adjusted to simplify For identical particles, the choice 
A = (2r3 — Vi — ra) / simplifies the treatment of permutation symmetry. 
For four particles, one can use 

miri + m2r2 ELi "^i^* /Ao^^ 

X(xr2-ri , y (xrs ■ , zocr^- , (A24) 

nil + ni2 nii 



© 




© 



FIG. 37: Jacobi variables for three particle systems. 




FIG. 38: Possible choices of Jacobi variables for four particles. 

or 

msrs + niir^ miri + m2r2 

X (X r2 — Vi , u (X — , v oc , ( A25) 

ms + m4 mi + 1712 

as illustrated in Fig. IHHl 

A more systematic way of constructing Jacobi coordinates consists first of rescaling = 
3^i/y?77^, and Pj = ^/rr^iq^, so that the kinetic energy reads {q1 + ■ ■ ■ + q1f)/2. Then, any 
orthogonal transformation leaving invariant the total momentum defines a possible choice 
of Jacobi coordinates. 

The choice of the Jacobi coordinates is often dictated by the approximation scheme one 
intend to use. For instance, if a four-body system is described as a weakly bound atom- 
atom system [(1,2) (3,4)], the choice ()A25j) seems convenient, with a wave function of the 
form = (p{x)(^{u)^(v). 

Note that for the Gaussian parametrisation described in Appendix El 

N 

^ = ^ exp[-qA('^)q/2] , qAq = , (A26) 

a change of Jacobi variables can be taken into account by a mere transformation of all 
matrices A^"-\ and thus should not modify the quality of the approximation. 



10. 2D binding 

It is well known that an attractive potential {V < 0, but V not always vanishing) in one 
dimension always supports a bound state with E < 0. For instance, one may consider V 
deeper than a tiny square well {V = for |x — Xo| > a, V = —Vq < if |x — XqI < a}, which 
is explicitly shown to support at least a bound state, however tiny are the size a and the 
depth Vq. 



The same result holds in two dimensions. See, e.g., Refs. |239l 240l |. A more precise 



statement is that V supports at least one bound state for n = 1 and n = 2 dimensions if 
/y(x)d"x < 0. 

A possible way of being convinced that n = 2 is the critical value of the dimension 
for binding in any attractive potential relies on the simple "delta-shell" interaction V = 
—gS{\x\ — R). There is a bound state if the zero-energy wave function has a node. 

One can choose the zero-energy wave function as 4>{x) = 1 for r = \x\ < R. At r = R, it 
receives an admixture of the other isotropic solution of A(f) = 0, and becomes, for r > R, 

<P{x) = 1 - g{r''" - R'~-) , (A27) 

with a logarithmic term in the limit where n = 2. For n > 2, r^~" has finite variation 
between r = R and r = oo, and thus requires a minimal strength g to direct (p toward the 
negative region. For d < 2, the negative region is reached, however weak the strength g. 



11. Potentials with asymptotic Coulomb attraction 

It is known that any Coulomb potential V = —g/r {g > 0) applied (we are back in 3D) to 
a particle of mass m gives an infinite number of bound states, with energy En = —mg'^/ (2n^), 
where n = 1,2, . . .. 

Besides r = and r = oo, the reduced radial wave function has n — 1 nodes located at 
r = Tj-"^ {i = 1, n — 1), and it is observed that rl^}i ^ oo as n ^ oo. 

Now the (n s) state can be considered as the ground state of a potential with V = —g/r 
for r > rfj"_\ and V = +oo at r < r^li- Hence a potential with an attractive Coulomb tail 
and a hard core of radius R always supports bound state, since one can arrange n such that 

e\ > R- 

Therefore, a potential which is anything below r < R, and exactly Coulomb beyond R 
will support bound states with deeper binding than in the case of a hard core of radius R. 

One can thus understand that a potential with a Coulomb tail and an arbitrary short- 
range part decreasing very fast beyond some radius R, will support bound states. 



more that 



APPENDIX B: VARIATIONAL METHODS 

The principle of a variational method is to write down an empirical expansion of the wave 
function on a set of basis functions involving free parameters and to adjust the parameters 
to optimise the results. 

In this section we overview the most popular forms of variational basis functions used 
in correlated few-body calculations and compare their relative efficiency and accuracy on 
selected few-electron systems. 

We also extend our discussion to hyperspherical harmonics, numerical grid, Monte-Carlo 
methods, and compare them to conventional variational methods in terms of convergence 
and accuracy. 

1. Parametrisation of the wave function 

With some appropriate form of the basis functions, appropriate (e.g., exponential, har- 
monic oscillator, Coulombic or Gaussian), the matrix elements can be calculated analytically. 
Expanding the wave function in such a basis leads to an a accuracy of the method that can be 
improved almost arbitrarily by increasing the number of basis functions. The convergence, 
however, depending on the form of the basis functions, is often very slow. 

a. Hylleraas and Exponential basis functions 

For small systems (two- or three-electron systems) traditionally the best results have 
been obtained by variational calculations in Hylleraas coordinates in which the trial wave 
function for the example of a two-electron system is written in the following form 

exp(-ari -/^ra) J] aijks'uH'' [ , (Bl) 

i,j,k ) 

where u = |ri — r2|, t = r2 — ri and s = ri+r2- A is the antisymmetriser, aijk are the linear, 
a and (5 are the nonlinear variational parameters. The convergence is controlled by VL. 

The results can be significantly improved by combining linearly two or three wave func- 
tions of the above type with different values of a and (5 j4^. The different values of a and 
j3 control the asymptotic, intermediate and short range sectors: 

{p i+j+k<n ^ 

^exp(-a„ri-/?„r2) Yl 41''^''^ (■ (^2) 
71=1 i,j,k ) 

Pekeris et al. |24l| improved the Hylleraas functions by including terms which are essential 
for the short-range behaviour ( (rf + r|)^/^ — » 0) 

vi/(ri, rs) = ^{exp(-«n - iSr^) a,,kin.s' uH\s' + e)"\\nsr} . (B3) 

i,j,k,l,m 



Thakkar and Koga j242| later introduced an other variant of the above ansatz with 
fractional powers. 

A different exponential basis, without the explicit inclusion of the powers in the radial 
coordinates, is also very powerful [24rii] 

i<N 

^'(ri, ra) = aiexp{-airi - /^^rg - 7*^12)} , (B4) 

i 

where a^, /5j and 7j are chosen in a quasirandom way, so the calculation is similar to a 
Monte-Carlo calculation with a random distribution of the exponential factors. 

Korobov significantly improved Thakkar's method by int rodu cing complex aj,/3j,7j 
parameters making the variational ansatz more flexible. Frolov j244l | used the similar vari- 
ational function but carefully optimised the nonlinear parameters instead using random 
distributions. 

These forms represents a different direction in the computations. In the early days when 
the computational power was very limited, compact analytical representations were essential. 
With the advance of computers simpler basis functions (but larger basis dimensions) became 
more popular. 

The exponential/Hylleraas type functions are very powerful but the analytical evaluation 
of the bas is fun ctions is complicated or even impossible for more than two or three electrons. 
Goldman |245j | suggested a quite a different approach: his idea is to uncouple the correlations 
by using appropriate new variables. To calculate the radial part of the matrix elements of 
the Hamiltonian one has to integrate 

/•OO /"CXD /'Tl-\-r2 

1= ndri / rsdrs / ri2/(ri, ra, ri2)dri2 , (B5) 

Jo Jo J\ri-r2\ 

where /(^i, r2, ri2) comes from the basis functions and Hamiltonian. By introducing the 
perimetric variables 

s = ri+r2, v = ri2/s, w = (ri - r2)/ri2 , (B6) 
the integral can be rewritten as 

00 />! nl 



sMs / v'dv / (1 -fW)/(s,t;,w)dw , (B7) 
Jo Jo 

and if the / function has the form 

/(s, V, w) = S{s)V{v)W{w) , (B8) 

then the integral is fully uncoupled into a products of one-dimensional integrals. With the 
help of this trick certain forms of correlated basis functions can be used for systems with 
larger number of electrons. 

Signiflcant progress has also been made in handli ng integr a ls wi th an exponential on the 
relative distances for the four-body problem. See j246l 247 . 248l |. and references therein. 
This offers an alternative to the Gaussian parametrisation. 



b. Gaussian and correlated Gaussian basis 



The Correlated Gaussian (CG) of the form 

1 1 ^ 

GA{r) = exp{--rylr} = exp{-- ^ij^i ' ^j} (B9) 

is a very popular choice of basis functions in atomic and molecular physics. Here f stands 
for a one-row vector whose ith element is r^. The merit of this basis is that the matrix 
elements are analytically available and unlike other trial functions (for example, Hylleraas- 
type functions) one can relatively easily extend the basis for the case of more than three 
particles. The well-known defects of this basis are that it docs not fulfil the cusp condition 
and its asymptotics does not follow the exponential falloff. This latter problem, especially 
for bound states, can be cured by taking hnear combinations of adequately chosen CGs. The 
translational invariance of the wave function is ensured by requiring that the parameters A 
fulfil some special conditions. These conditions ensure that the motion of the centre-of-mass 
is factorised in a product form. 

By combining the CG with the spin function parts, the full basis function takes the form 

^kLS^A{xsMsGA{r)}, (BIO) 

with an appropriate spin function Xsmsj where "A;" is the index of the basis states and A is 
an antisymmetriser for the identical fermions. 

Instead of optimizing the parameters of A it is more advantageous to rewrite the equation 

as 

exp{ - ^ - r^)' - (Bll) 

i<j i 

The relationship between aij,Pi and A is 

aij = -Aij {iT^j), (3i = Y,^ki, (B12) 

k 

where aji {i < j) is assumed to be equal to aij. There are two reasons to choose this form. 
The first is that in choosing aij in this way we deal with a correlation function between the 

particles i and j, while A^j has no such direct meaning and during the optimisation it is 
more difficult to limit the numerical interval of A^j to be chosen from. Secondly, one can 
utilize this specific form to make the individual steps of the parameter selection very fast. 

2. Optimisation of the parameters 

To obtain very precise energy, one has to optimize the variational parameters Uki and 
Akij of the trial function. The dimension of basis sets is typically between 100 and 1000, 
and each basis state has nine nonlinear parameters. The optimisation of a function with a 
few thousands nonlinear parameters cannot be done efficiently by using a straightforward 
deterministic optimisation method, since this could entail the complete reconstruction of the 
Hamiltonian matrix and diagonalisation every time when some of the nonlinear parameters 
are altered. Moreover, the deterministic search for the optimal value of such a large number 
of parameters is likely to get trapped in a local minimum. 



Different strategies have been proposed. Examples are given below. A comparison of 
results obtained for the celebrated benchmark case of the He atom is proposed in Table . 



TABLE V: Variational ground state energy of the He atom obtained by various trial functions. 



Method Ref. Basis size 


Energy (a.u.) 


Pekeris [241] 


246 


-2.9037243770340 


Thakkar [242] 


308 


-2.9037243770341144 


Drake [24^] 


1262 


-2.90372437703411948 


Goldman [245] 


8066 


-2.903724377034119594 


Korobov [33] 


2114 


-2.903724377034119598296 


Drake [4^ 


2358 


-2.903724377034119598305 


Baye M 


18900 


-2.9037243770340 



a. Frolov's method 

This author mainly used exponential functions to perform a variety of three-body calcu- 
lations and some studies of larger systems. His so-called "two-stage" strategy is presented, 
e.g., in Ref. ^25 L] . Only the first non-linear parameters are optimised exactly. The next 
terms in the variational expansion are chosen to vary in a regular manner that is described 
by just a few parameters. 



b. Kamimura's method 

The group of Kamimura has used Gaussian wave functions in a variety of few- body 
calculation. A comprehensive description of the method is given in a review article |252j l . 
There are two main ideas. Firstly, in a given variable r, the range parameters ai of ip oc 

Ci exp(— Q;jr^/2), instead of been independent, are assigned to belong to a geometric series, 
say aj+i/ttj = = ■ ■ ■ . Hence, only the first and last terms have to fitted. Secondly, 

for more than two Jacobi variables, say x and it is not attempted to write down the most 
general Gaussian exp[— (aa;^ + 2f3x.y + 7y^)/2]. Instead, if x describes the relative motion 
inside the (1,2) cluster, and y the motion of the third particle relative to this cluster, one 
seeks at a component 

c^ exp[-a,ajV2] \ [Y.'^J exp[-/3,t/V2] J , (B13) 

in the wave function, to be supplemented by analogous terms with clusters (1,3) and (1,2). 
The physical interpretation is to construct interferences of several virtual-decay channels. 



c. Stochastic variational method 



A procedure based on the stochastic search for the best set of nonlinear parameters can 
be programmed efficiently and is capable of achieving highly accurate results for most few- 



body systems. The essence of the strategy can be summarised as follows: Let {ui,Ai} be 
the nonlinear parameters of the ith basis function out of the set of K such basis functions. 
Then the procedure is 

(1) A succession of different sets of {{u},A}}, {u^' , A^'}) are generated randomly. 

(2) By solving the eigenvalue problem, the corresponding energies {Ej , E^") are deter- 
mined. 

(3) The parameter set {u^, A^} which produces the lowest energy is then used to replace 
the existing {ui, Ai} set. 

(4) The procedure cycles through the different parameter sets {{ui,Ai},i = 1,...,K), suc- 
cessively choosing different sets to minimize the energy until convergence is reached. 

The essential reason motivating this strategy is the need to sample different sets of nonlinear 
parameters as fast as possible. The main advantage is that it is not necessary to recompute 
the complete Hamiltonian nor it is necessary to solve the generalised eigenvalue problem 
from scratch each time a new parameter set is generated. By changing the elements of 
parameter set for each basis function individually, it is necessary to recompute only one row 
(column) of the Hamiltonian and overlap matrices each time the parameter set {ui, Ai} is 
changed. Furthermore, the solution of the generalised eigenvalue problem is also expedited 
since the Hamiltonian matrix is already diagonal apart from one row and one column. 

A similar strategy to the above was used when adding additional terms to the basis. 

The above way of finding the best parameters is certainly very restricted. E ven t his simpl e 
method gives very accurate energies, as seen in applications to various fields jlslBliH^ . 
More sophisticated techniques may give better results in a smaller basis size. 



3. Real space approaches 

The representation of wave functions in "real space" on a numerical grid (mesh) is a very 
simple and powerful technique. In the simplest case, the basis functions are simple Dirac 
delta functions 

Mx,) = 5,, (B14) 

in the coordinate space point Xj. A more sophisticated choice pioneered by Baye et al. j25nj | 
uses infinitely differentiable Lagrange-Laguerre functions 

/.(x,) = {-irx\"^^^eM-x,m (B15) 

j i 

where Ln{x) is the Laguerre polynomial of degree and Xi is one of its zeros, i.e., L^{xi) = 
0. The fi{x)fj(y)fk(z) product of the one dimensional basis functions can be used to expand 
the functions of relative motions in an A^-electron problem. This representation is very 
flexible and approximates the exact wave functions extremely well. 

The advantage of the mesh techniques is that matrix elements can be trivially calculated. 
The potential matrix is diagonal therefore it is particularly simple independently of the form 
of the potential. The Hamiltonian matrix is very sparse so powerful linear algebraic iteration 
techniques such as Conjugate Gradient method with efficient multigrid preconditioning can 
be used to solve the eigenvalue problems. 



4. Hyperspherical harmonics approach 



The method of hyperspherical expansion is wi dely used in atom i c phy sics, especially in 
its variant "adiabatic approximation". See, e.g., fl^L I255L l25d 257 . l258| |. and the detailed 
review by Lin j2^ for details and references. In the accurate calculation by Mandelzweig 
et al., the method is compared with other variational methods (t^ . 

Schematically, in the case of three particles, the two Jacobi variables p and A are replaced 
by an overall distance and a set of five angles, 

r = (p2 + , ^^5 = {p, A, tan-i(p/A)} . (B16) 

A partial wave-function expansion 

[L] 

in terms of the generalised harmonics y results into a set of coupled differential equations 
for the radial components Here, [L] denotes the "grand" angular momentum L and 
associated magnetic numbers compatible with rotation and permutation invariance. The 
convergence is studied as a function of the maximal value of L which is allowed. In practice, 
the number of equations grows very rapidly as L is incr eased . Several approximation schemes 
have been proposed, such as "potential harmon ics" \25l\ and "adiabatic approximation" 
255l 259^ . An interesting variant is proposed in ||26C| . 



5. Faddeev equations 

The Faddeev equations |26l| , first written to describe scattering in momentum space with 
short-range potentials, have been reformulated in position space, and applied to a varietv o f 
problems involving even confining potentials or long-range Coulomb potentials |262L 1263^ , 
after suitable modifications. 

If the Hamiltonian reads 

H = T + Vi + V2 + V3 , (B18) 
with the kinetic energy T and the pair potentials Vi = f (r23), etc., then a decomposition 

^ = ^1 + ^2 + ^3, (B19) 

associated to the set of equations 

{E~T)^i = Vi^, z = l,2,3, (B20) 

implies the original Schrodinger equation H'^ = E^^. Permutation symmetry is easily 
implemented by imposing relations among the 

In practice, one often used a partial- wave expansion of each Faddeev component \E'j(p, A) 
to account for all internal orbital momenta ip and i\ compatible with a given overall angular 
momentum J. One obtains a set of coupled integro-differential equations in p and A. The 
convergence is studied by varying the maximal value allowed for ip and ix. 

Faddeev equations are not primarily aimed at competing to obtain the most accurate 
values of the binding energy. Their main merit is to produce high-quality wave-functions, 
accounting for the correlations, even when restricted to the lowest values of ip and ix. 



6. The Variational Monte Carlo method 



The Variational Monte Carlo method (VMC) is a very powerful numerical technique that 
estimates the energy and all desired properties of a given trial wavefunction without any 
need to analytically compute the matrix elements. It poses no restriction on the functional 
form of the trial wavefunction, requiring only the evaluation of the wavefunction value, its 
gradient and its Laplacian, and these are easily computed. Using the VMC algorithm, 
essentially a stochastic numerical integration scheme, the expectation value of the energy 
for any form of the trial function can be estimated by averaging the local energy H^> /^! over 
an ensembl e of configuration s dist ributed as \E'^, sampled during a random walk using the 
Metropolis [iei or Langevin |265j algorithms. The fluctuations of the local energy depends 
on the quality of the trial functions, and the VMC can also be used to optimize the trial 
functions. 

A popular and effective approach to building compact explicitly correlated wavefunctions 
is to multiply a determinant wavefunction by a correlation factor, the most commonly used 
being the Jastrow factor. The inclusion of the Jastrow factor is only possible by VMC 
techniques. 

There are various sophisticated trial wave functions have been implemented and tested. 
The most simple but typical trial function is written assuming a Pade factor exp[(ar + 
6r^)/(l + cr)] for the electr on nu cleus wave function and a Jastrow factor exp[a'r/(l + c'r)] 
for the interelectronic part l2f 



E 



arj + hr] 

1 + CTi 



exp 



e 



SM 



(B21) 



where Qsm is the spin function. The Pade factor is simple but good choice for the nucleus- 
electron part providing flexibility with small number of parameters, representing the coales- 
cence as well as the decay of the wave function for r — and r ^ oo, respectively. 

As an example, we compare in Table IVll the ground state energy of the Be atom obtained 
with this simple VMC ansatz with results of other approaches including those which use 
correlated Gaussian basis functions. This simple form does not give very precise energy, 
but the energy can be improved by using more elaborate correlation factors (an example is 
shown in the second line of Table rvl|l or by using diffusion Monte-Carlo method. The main 
advantage of the VMC is that while it is not necessarily the most precise method for small 
systems its applicable to larger systems with favorable scaling properties. 



TABLE VI: Variational ground-state energy of the Be atom obtained by various trial functions. 



Method Ref. Basis size Energy (a.u.) 
VMC 1 \2M -14.6528 
VMC 2 [266] -14.6651 
CG [267] 1200 -14.6673 
SVM [254] 500 -14.6673 



7. Short-range correlations 



A quantity of interest for exotic systems involving positrons and antiprotons is the lifetime 
due to internal annihilation of a (e"*", e~) or a (p, p) pair. It is proportional to the probability 
per unit volume of finding these constituents at zero separation, i.e., to an expectation value 
of the type 

6., = {mr^-r,)\^) , (B22) 

within the exact wave function assumed here to be normalised. 

Computing 6ij accurately is much more difficult than estimating the binding e nergy , 
or an observable that is spread over the whole range of the wave function It^ . l268l | . 
It often happens that an approximate variational wave function is excellent in the region 
relevant for estimating the energy E = (H), but, for instance, underestimates the wave 
function at large distances, and then overestimates its value at large separation, in order 
to remain properly normalised. This is the case when harmonic oscillator wave functions 
^ c„if„(Q;r) exp(— a^r^), where Hn is a polynomial, are used to describe binding by non- 
confining forces. The Gaussian parametrisation ^ 7^ exp(— /^jr^) is better in this respect, 
thanks to the freedom of introducing a variety of range parameters. 

Note that if a sequence of approximate normalised wave functions \E'j leads to energies 
Ej that converge toward the exact energy E, it is only guaranteed that, if $ denotes the 
exact normalised wave function, 

(B23) 

but not that \E'j $ (modulo a phase) at any point. 

To obtain an accurate estimate of 6ij without pushing too far t he variational expansion 
for the one can use a trick devised by Schwinger (see, e.g., |238l|), and further developed 
for estimating parity-violating effects in atoms. The basic idea is that if 

= ($151$) , (B24) 

for the exact wave function, then the sequence 

B, = {^iI,\B\^,) , (B25) 

of approximate values might converge much better than the plain Aj = {'^jlAl'^j) to estimate 
{A). 

For N = 2 constituents, if the wave function is written as usual $ = Yf^{Q)u{r)/r, with 



u"ir) - r^l±llM(r) + - \/(r)]M(r) = , (B26) 



the Schwinger rule reads |238l | 

= A r[V\r) - 2i{i+l)/r']dr , (B27) 
4vr Jo 

and turns out very powerful. In particular with a linear potential sometimes used for de- 
scribing quark confinement, the right-hand side always gives the right answer for i = 0, 
however poor is the variational approximation for u{r), provided it is normalised. For i > 0, 
612 vanishes, and the rule ()B27|) links the expectation values of V'{r) to that of r~^. In the 
Coulomb case, there is, indeed, a well-known identity between (r~^) and (r~^). 



The generalisation of the Schwinger ru le to systems with > 2 constituents has been 
written down by Hofmann-Ostenof et al. j269l | and by Hiller et al. j270j . For constituents 
with mass m = 1, it reads 



An 



2£? 



where li = xixp^ is the angular momentum carried by the first Jacobi variable Xi 
and the differentiation Vi should be understood with the other Jacobi variables X2, 
kept fixed. 



(B28) 
r2-ri, 



Evaluating the right-hand side to improve the mere reading of 6ij from variational wave 



functions, was us ed t o study parity- violating effects in atoms ^70j, or short-range correla- 
tions in baryons 262l |. Consider for instance Ps2. A very crude approximation of the wave 
function of Ps2 gi ves a probability 5(e"'"e~) = 1.84 x 10^^ by direct reading, and is corrected 
into 2.16 X 10~^ |l83l | by the g enera lised Schwinger rule, close to the value 2.19 x 10~^ 
obtained by refined calculations |l5C| . 

Several variants for calculating 5ij have been proposed |27lj . Perhaps, the method of 
matrix-element identities works best at the beginning of a variational expansion, to speed 
up the convergence of 6ij towards the neighbourhood of the ex act v alue, but does not help 
much once one is seeking at very high accuracy. According to |272| . indeed, "for low-lying 
states, there is no substitute for direct high-precision calculations of the wave-function near 
the origin." 



APPENDIX C: THE BORN OPPENHEIMER APPROXIMATION 



In 1927, soon after the discovery of quantum mechanics, Born and Oppenheimer |273l | 
were able to obtain an approximate solution of the Schrodinger equation for a molecule by 
separating the electronic and nuclear motions. This was possible on account of the much 
larger masses of the nuclei, as compared to the mass of the electron. Born and Oppenheimer's 
treatment can be formally justified in terms of a perturbation expansion of the exact non- 
relativistic energy and wave function in terms of the parameter, (me/M)^/"^, where mg is the 
mass of the electron and M is the average mass of the nuclei. 

This approximate quantum mechanical treatment of molecules, known as the Born- 
Oppenheimer (BO) approximation, has always played a central role in the study of molecular 
structure. Descriptions of the BO app roxim ation are given in, for example, Bransden and 
Joachain ^| and Pauling and Wilson |l44| . 

For simplicity, let us consider how the BO approximation is applied to a diatomic 
molecule, containing N electrons and two nuclei, A and 5, with masses Ma and Mb and 
charges Za and Zb-, respectively. Our treatment can be extended to more complicated 
molecules. 

The first step in the BO approximation is to fix the nuclei in space with a given R value 
and solve the Schrodinger equation for the electrons interacting with each other and with 
the nuclei through the Coulombic force. The equation is of the form 

H,,^{r-R)=E{R)m{r-R) , (CI) 

where He\ is the electronic Hamiltonian. \l/(r; i?), where r = (ri, r2, . . . , tat), is the elec- 
tronic wave function and E{R) is the sum of the electronic energy at the given R value and 
the energy due to the Coulombic repulsion of the nuclei. 
Hei is of the form 

N N N N 

^ i=i i,j=i 1=1 i=i 

where is the distance between electrons i and j and and ViB are the distances of 
electron i from nuclei A and B, respectively. 

The electronic wave function, \l/(r; R), depends parametrically on the internuclear vector, 
-R. Its dependence on the direction of R is easily taken into account. The z-axis of the 
electronic coordinates is normally taken to be along R; the probability of a given orientation 
is determined at a later stage with the calculation of the nuclear wave function, x{R)- 

E[R) and \l/(r, R) are determined as functions of R. Except in very simple cases such as 
H^, for which N = 1, for which equation ()C1|) is separable in prolate spheroidal coordinates, 
only approximate solutions can be obtained using the Rayleigh-Ritz variational method. 
The quality of the approximate solutions for a given molecule has improved very significantly 
over the years with increases in computer power and the ingenuity of quantum chemists. 

The nuclei are considered as two particles moving in the central potential, V{R), where 

V{R) = E{R) . (C3) 

Thus their centre of mass motion is separated out in the usual way. The nuclear wave func- 
tion, x(-R); for the internal motion of the nuclei, is determined by the nuclear Schrodinger 



wave function, 

HnX{R) = EtX{R) (C4) 

where 

Hn = -^Vl+ViR) , (C5) 

and 

MaMb ,„„, 

is the reduced mass of the nuclei. Et is the energy, in the BO approximation, of the molecule 
under consideration. The associated BO wave function, \E'T(r,i?) is of the form 

^T{r,R) = x{R)'^{r;R) . (C7) 

The most common states to be treated using the BO approximation are the states as- 
sociated with the electronic ground state, i.e. the solution, (or solutions, in the event of a 
degeneracy), of equation ()C1|) that corresponds to the lowest minimum value of E{R), as 
a function of R. However, excited electronic states can also be treated. In principal, it is 
possible to obtain a complete set of BO solutions of the form in equation ()(^7j) . This would 
necessarily involve continuum states. 

In the BO approximation the terms in the full Schrodinger equation for the internal 
motion of the molecule that represent coupling between the nuclear and electronic motion 
are neglected. They are the term containing 

Vnx{R)-'^R'^{r;R) (C8) 

and the term containing 

x{R)V],^{r;R), (C9) 

which result from the operation of Vi? on the electronic wave function. 

Also the ratio of the mass of the electrons to the masses of the nuclei is assumed to be 
infinite, resulting in their neglect in the centre of mass separation. This would have to be 
corrected in an exact treatment. Allowance for the large but not infinite electronic- nuclear 
mass ratio leads to the introduction of small mass polarisation terms in the Hamiltonian that 
couple the momenta of the electrons and also, if 7^ Mr, the momenta of the electrons 
and the nuclei. See, for example, Kolos and Wolniewicz 2741 . 



A simple improvement to the BO approximation, called the adiabatic approximation, is 

to set 

V{R) = E{R) + A{R) , (CIO) 

where 

A{R) = (vl/(r; R) \ -i- | ^{r; R)) . (Cll) 

In equation ()(^11|1 . the round bracke ts iii dicate integration with respect to the electronic 
coordinates. See, for example, Kolos j274l | and Messiah |13]. Further improvements to the 
BO approximation are considered by Pack 275| |. 



Calculations of corrections to the BO approximation as applied to and IID+ are 
referenced in Sec. IIII Bl The corrections to the BO approximation for the groui id st ate of 
diatomic two-electron molecules have been calculated by Kolos and Wolniewicz EtI . The 



correction to the energy is small in comparison with the binding energy of these molecules. 



4.74 eV, calculated using the BO approximation. This is a consequence of the large value 
of Ai/me > 918.1. 

However, as pointed out in Sec. IIII Gl the BO approximation does not work well for 
muonic molecular ions. In the case of dt/i, for example, the ratio of the reduced mass of the 
d and t to the mass of the electron is 10.7, which is much smaller than the value of /J, /me 
in the case of diatomic, two-electron molecules. 

A co mpa rison is made of the binding energies of the five bound states of dtfi in the table 
in ref. |276l | as calculated using the BO approximation, the adiabatic approximation and 
by a very accurate variational calculation using the full non-relativistic Hamiltonian for the 
internal motion of dtfi. The binding energies for the most weakly bound state, the J = 1, 
V = 1 state that plays a crucial role in muon catalysed fusion, obtained by the BO and 
adiabatic approximations and the variational method are 9.7, 7.7 and 0.66 eV respectively, 
showing the considerable errors that result from the two approximations. 

It is of interest to note that the adiabatic represent ation method referred to in Sec. IIII G\ 
which was used by Ponomarev and his co-workers |l09j | to correct the error in the BO 
approximation, involved the use of many BO basis functions of the form in equation ()(y7|l in 
a close-coupling type calculation using the full Hamiltonian for dt/i referred to above. For 
details of applications of refinements of the BO a ppro ximation to the t/i -|- D2 reaction in 
the muon catalysed fusion cycle, see Zeman et al. [277|. 

Co rrect ions to the BO approximation in the case of H H ha ve been calculated by Armour 
et al. |213j using the method described by Hunter et al. 

A very similar approximation to the BO approximation, the adiabatic-nuclei or nuclear- 
impulse approximation, can be used to separate the nuclear and electronic mo tion in scat- 
tering calcul ation s invol ving atoms and/or molecules. See, for example. Chase |279| | , Ch ang 



and Temkin [280[, Lane [281|. For applications to positron scattering, see Armour 



APPENDIX D: ELEMENTARY THREE- AND FOUR-BODY CALCULATIONS 



1. Introduction 

In this Appendix, we present a few basic calculations concerning systems of 3 or 4 charges. 
Some of them are borrowed from pioneering papers that illustrate the ability of quantum 
theory and variational methods to reproduce or even predict the properties of systems more 
complicated than the hydrogen atom. The analytic results presented below, when properly 
generalised, are basic ingredients of powerful variational methods presented in Appendix O 



2. Helium and Helium-like atoms with fixed nucleus 

This corresponds to the Hamiltonian 

H=f + ^ + — , Dl 

where the electron mass is set to m = 1. A rescaling is often used to set the attraction 
strength to unity and the repulsion strength to 1/Z, leading to the well-known ex- 
pansion. See, for instance, 136l 137l |. 



When the repulsion among electrons is neglected, one can solve exactly the unperturbed 
Hamiltonian Hq = H — l/ri2 and obtain for the ground state an energy Eq = —Z^ and a 
wave function \E'o oc exp(— Zri — Zr2). 

Assuming ^ to be normalised, the repulsion gives a first order correction 

Eo-^Eo + {^o\r^2\'^o) . (D2) 

The integral corresponding to the repulsive term is often cal culat ed by means of an expansion 
of r^2 terms of spherical harmonics. Se e, for instance, 283l |. For the ground-state, one 



can also use the Gauss theorem, as, e.g., in |284| . The field created at distance r by the first 
electron is 

1 1/1 2Z \ 

^(^) = — 2 ^(^)' = — ^ + ^ + — + exp(-2Zr) , (D3) 







^2 \ ry'2i 



where u{r) = 2Z^/^r exp(— Zr), corresponding to a potential 

V{r)= E{x)dx = — / u{xfdx- dx 

Jr Jo J r ^ 

= ^ +(^ + z\ exp(-2Zr) , 



(D4) 



and thus a repulsion energy 



{ru') = -l V{r)u{rydr = ^. (D5) 



To first order, the energy of helium is thus —11/4 = —2.75. 



A standard improvement consists of using a normalised wave function 



oc exp(— ari — ar2) , (D6) 

with kinetic en ergy and potential energy —2Za for the attractive part, and 5a/8 for the 
repulsive part |285l |. The parameter a is interpreted as the effective charge seen by each 
electron. One should thus minimise the variational energy 

E{a) = c? - 2Za + 5a/8 = [a-{Z - 5/16)]^ - {Z - 5/16)^ , (D7) 

this giving an improved minimum E = —729/256 ~ —2.85 at a = 27/16 ~ 1.69, in the case 
of Z = 2. 

Note that the threshold for dissociation of charges {Z, —1, —1) into an ion {Z, —1) and 
an isolated electron is Eth = —Z'^/2, so that miiiaE — Eth = —{Z — 5/16)^ + Z^/2 remains 
negative only for Z > (10 + 5\/2)/16 ~ 1.07. For smaller values of Z, demonstrating the 
stability of [Z, —1, —1) requires trial wave functions more elaborated than ()D6|) . This is 
in particular the case for H~, where Z = 1. The stability of H~ cannot be reached using 
factorised wave function /(ri)/(r2), however refined the individual function /. 

Looking more closely at H~, again in the limit where t he proto n mass is infinite, the 
simplest wave function which achieves binding is [2(1 licl li^ l6l[ Il30| |. 

*[a, b] = exp(-ari - fora) + {1^2}, (D8) 

or a factorised wave-function with a = b but an explicit account of the correlation, e.g., by 
a factor (1 + cr^), or a combination of a b and of the ri2 term. For the choice ()D8|1 . the 
matrix elements of normalisation, potential energy and kinetic energy are easily calculated 

27r2 27r2 



2, 7r2(a2 + 62 2a67r2 



^^^^ ^ cfi ' (D9) 

(r-i) =7r2(a + 6)(a-=^r=^ + a-6), 
, 7r2(a2 + 3a6 + 62) ^t^^ 



where a = {a + b)/2. 

In the case of Helium, one obtains an energy E = —2.87566, i.e., a modest gain with 
respect to the case where a = b {E = —2.8476). However, the values a ~ 1.19 and b ~ 2.18 
of the parameters show a non-negligible anticorrelation between the two electrons. 

In the case of H~, binding is reached since the variational energy E = —0.513303 lies 
below the threshold -Eth = —1/2 by a fraction x =~ 0.027. Anticorrelation is very strong, 
as 6/a ~ 3.67. This means that if one electron is close to the proton, the second is far away, 
and vice-versa. 

The Hughes-Eckart or mass-polarisation termp^.p2 has obviously a vanishing expectation 
value within this wave function ()D8|) . Thus, if the finite mass of the nucleus is taken into 
account, both the threshold energy and the variational energy are rescaled by the same factor 
M/{M + m), where m is the mass of the electron and M that of the nucleus. Thus, with 
this wave function, the fraction of binding is the same for all isotopes of Helium (neglecting 



finite-size effects) as for a static Z = 2 source. Similarly, as noticed by Hill |20|, all systems 
{M~^m~m~) are demonstrated to be stable (and within this approximation, with the same 
fraction x and same anticorrelation b/a). 

The wave function ()D8|1 or even a superposition of wave functions of this type will never 
approach very closely the exact solution, as it lacks one degree of freedom, namely some 
dependence on ru, the distance among the electrons. A wave function that allows for this 
ri2 dependence is 

\l/[a, b, c] = ip[a, b,c] + ■ ■ ■ , ip = exp{—ax — by — cz), (DIO) 

where the dots denote terms deduced by symmetry. Here x = r23, etc. For the matrix 
elements deahng with scalar states, the integration is over xdx ydy zdz, subject to triangular 
inequalities. A basic integral is 

F[a, b,c] = J exp(— 2ax — 2by — 2cz)dxdydz, 

oo oo +x 

= \ j exp(-2ax)da; j exp[-(6 + c)a]d(j j exp[-(6 - c)5]d5, (Dll) 



4(a + b){b + c)(c + a) 
in term of which the normalisation of a single exponential ip is 



and a typical term of the potential energy is 

1 dF[a, b, c] 

For those local operators, off-diagonal matrix elements between ip[a, b, c] and ip[a', b', c'] are 
deduced by the substitution a —>■ a = {a + a')/2, etc. 

For the kinetic energy, p^ip = —iip{by — cz), and thus the most general matrix element is 

{ip[a',b',c']\pl\(l)[a,b,c]) = {bb' + cc')n{a,b,c) + {be' + b'c)ti{a,b,c), (D14) 

where ti is the matrix element of 

ii = y.z = "^'-f-'" ^ (D15) 
2yz 

and is given by a combination of suitable derivatives of the generating function F[a, b, c]. 

For the Helium atom, the wave function ()D10|) contains two exponentials related by 
{1 2} exchange and thus three parameters a, b, and c (two, if one use the virial theorem). 
It provides a variational e nerg y E = —2.89953, to be compared with the best available result 
E = —2.90372 (see, e.g., |259l | . and references therein). 

For H~, the fraction of binding reaches Xb = 0.047696, while for Ps~, it is x = 0.0268. 



3. Four-body systems 

Let us consider here a system (1, 2, 3, 4) with charges ±(+l,+l,— 1,— 1). We shall only 
consider the groun d st ate. A wave function that is a generalisation of the one used by 
Hylleraas and Ore [l49| for Ps2 reads 

'^[a,b,c,d]=ip[a,b,c,d] + ^^^^^ 
(p[a, b, c, d] =exp[-(ari3 + br 23 + cru + dr2i)/2], 

where once again the dots means exponentials deduced from symmetry considerations. In 
matrix elements not involving r34, the integration runs over 

dr = ri3ri4r23r24dr', dr' = dri2dri3dri4dr23dr24, (D17) 

as the system is fully identified by the two independent triangles (1,2,3) and (1,2,4). A basic 
integral is 

[ dr' 

F[a, 6, c, d, ti] = / exp(— ari3 — 6r23 — cr 14 — (ir 24 — Mri2). (D18) 

J ri2 

Again, the triangle inequalities are conveniently accounted for by using the variables 
(c^s, S3) = ^"13 ± ^23, and similarly for particle 4. One obtains 

Fla, i, c. d, u] = " + " + ""t + " + f "° + ° + + + "»l . (D19) 

^ ^ {a-b){a + b){c-d){c + d) ^ ' 

The normalisation integral, and the matrix elements of the potential can be written as 
derivatives of F . For instance the normalisation is 



n{a,b^c,d)= I \(p^\dT = — ^ 



dudadbdcdd 



(D20) 



u=0 



The internuclear repulsion is 



_i\ 9^Fi2[a,6, c, d, 0] 



""'"■"■'■^l O^OWM ■ (°^') 
The electronic repulsion is obviously 

f 34 [a, 6, c, d] = f 12 [a, c, b, d] . (D22) 

The kinetic energy of particle 3 can be calculated starting from 

P3ip[a, b, c, d] = -iip[a, b, c, d]{af 13 + 6f23), (D23) 

so that, exactly as in the 3-body case 

{^\pI\^) = {a' + b'){v\v) - 2ab{^\{rl, - r^, - r223)/(2ri3r23)|^) (D24) 

which can be expressed in terms of derivatives of F[a,b,c,d,u]. The expectation values of 
pf for other particles are deduced by suitable permutations of a, b, c and d. 

Off-diagonal matrix elements between ip[a, b, c, d] and (f[a', b', c', d'] are given by the simple 
rule a — i> a = (a + a')/2 for the local operators. For the kinetic energy terms in Eq. ()D24|) . 
one replaces + 6^ by aa' + bb' and 2ab by ab' + a'b, while the operators adopt the average 
arguments a, b, etc. 

One eventually obtains explicit expressions for the expectation value of the normalisation, 
kinetic energy and potential energy. Equation follows straightforwardly. 
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